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SYNOPSIS 


TRIANGULAR SEGMENTS AND TWO-DIMENSIONAL ANALYSIS 
FOR MICROWAVE INTEGRATED CIRCUITS 
A Thesis Submitted 
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for the Degree of 
DOCTOR OF PHILOSOPHY 
by 

Rakesh Chadha 

Department of Electrical Engineering 
Indian Institute of Technology, Kanpur 
March, 1981 


The concept of two-dimensional planar circuits is a 
generalization of one-dimensional stripline or micro strip circuit! 
to the case when the transverse dimensions are comparable to wave- 
length. Thus, these circuits have dimensions comparable to (or 
greater than) the wavelength in two directions but have the third 
dimension much smaller than the wavelength. Analysis of two- 
dimensional circuits containing rectangular and circular segments 
by using impedance Green’s functions is available in literature. 
Investigations reported in this thesis are aimed at (i) analysis 
of triangular planar segments and their applications in microwave 
integrated circuits, and (ii) use of two-dimensional analysis for 
the design and optimization of strip line and microstrip circuits. 

The thesis begins with the development of impedance Green’s 
functions for the analysis of triangular segments. Three types of 
triangles (30° - 60° right-angled, equilateral, and right-angled 
isosceles) are considered by placing additional image sources out- 
side the triangular regions. The Green's functions satisfy 



magnetic wall condition all along the periphery of the segments. 
Green’s functions for some triangular shapes with electric wall 
on one of the sides and magnetic walls on the other two sides are 
also obtained. 

Green's functions for circular sectors, annular rings, and 
annular sectors have been obtained by expanding these in series 
of eigenfunctions. As mutually orthogonal eigenfunctions are 
obtainable only for those sectoral angles which are submultiples 
of 180°, the Green's functions are obtainable for only these 
sectoral shapes. 

Segmentation method for the analysis of planar circuits 
involves breaking up of a circuit pattern into segments for which 
Green's functions are known. S-matrix of the overall circuit is 
obtained from the S-matrices of the individual segments. This 
method is modified so that Z-matrices for individual planar seg- 
ments are used to obtain the overall Z-matrix. The proposed 
method using Z-matrices is shown to be computationally more effi- 
cient than the earlier method. 

The Green's functions, for the equilateral triangular seg- 
ments, derived in this thesis have been used to analyze rhombic 
and hexagonal resonant shapes by segmenting these into triangular 
shapes. Equilateral triangular, 60°-rhombic, and regular hexago- 
nal shapes are analyzed and their resonances have been investiga- 
ted. It has been pointed out that equilateral triangular shapes 
have odd modes of resonances in addition to the even modes 



w i 

reported earlier. For rhombic and hexagonal shapes, their reso- 
nance frequencies as well as the voltage variations around the 
periphery have been obtained for fundamental and some higher 
order modes. The resonance frequencies predicted have been veri- 
fied experimentally for each of three shapes. 

Stripline and microstrip T-junctions have been shown to 
exhibit lower parasitic reactances when triangular geometries are 
used as compared to the junctions with rectangular geometry nor- 
mally used. Two types of T-junctions with impedance ratios 1:1:1 
and ( 1 / V*2 ) : 1 : 1 have been studied for this purpose. For T- 
junctions with impedance ratio 1:1:1, an equilateral triangle is 
used at the junction. For T-junctions with impedance ratio 
(l/V2):1:1, a right-angled isosceles triangle (in magnetic wall 
model) is used at the junction. Improvements using triangular 
shaped junctions have been verified experimentally for both the 
cases . 

The two-dimensional analysis technique can be used in opti- 
mization of the design of the circuits whose performance is limi- 
ted by the junction reactances and the presence of higher order 
modes. This procedure has been used to obtain improved designs 
for 3 db in-phase power dividers and 3 db branch-line hybrids. 

Three types of power divider circuits that have been con- 
sidered are : i) a circuit with single matching section (Z=Z Q /f2) 
on the input side, ii) a circuit with single matching sections 
(Z=Z Q ][ 2 ) on the output sides, and iii) a circuit with a matching 



section (Z=Z 0 /2^^) at the input side and matching sections 
(Z=Z Q 2^^) on the output sides.. Triangular geometries have 
been used at the T-junctions . The lengths of various matching 
sections and position and value of isolation resistance are 
optimized for each of the three types of power dividers to 
minimize input and output VSWR's at the center frequency. 
Experimental results for a microstrip circuit with single 
matching section at the input side are included* 

For 3 db branch line hybrid also, right-angled isosceles 
triangles are used at each of the T-junctions, lengths of the 
quarter wave sections are optimized to minimize the input VSWR 
and to equalize the power division at the center frequency. 

These two circuit examples demonstrate the applicability 
of two-dimensional circuit analysis for the design and optimi- 
zation of strip line and microstrip circuits. 



Chapter One 


INTRODUCTION 

1.1 TWO-DIMENSIONAL COMPONENTS 

Two-dimensional (2-d) microwave planar components have 
been proposed for use in microwave integrated circuits 
[l] - [4] . These components have dimensions comparable to 
(or greater than) the wavelength in two directions but have 
the third dimension much smaller than the wavelength. Three 
types of configurations for the 2-d components are thus 
possible. These are -(i) triplate or stripline type shown in 
Pig. 1.1(a) and (c),* (ii) open or microstrip type shown in 
Pig. 1.1(b) and (c),* and (iii) waveguide or cavity type shown 
in Pig. 1.1(d). The triplate and the open type configurations 
can be considered to be generalizations of one-dimensional 
stripline and micro strip circuits to the cases when the trans- 
verse dimensions are comparable to the wavelength. The wave- 
guide type configuration (Pig. 1.1(d)) is a special case of 
three-dimensional waveguide in which the height is much smaller 
than the wavelength. 

1.1.1 Me thods of Analysis 

Various methods for obtaining the characterization of a 
2-d component depend upon its geometry (i.e. shape of the 
central conducting patch in a triplate structure or the upper 
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conductor patch in microstrip type circuit). When a 2-d compo- 
nent is of a simple geometrical shape, impedance Green' s fun- 
ctions [l] are used for its analysis. These Green's functions 
have been known for only rectangular and circular shapes 
[3]> [5]. Using Green's function, impedance matrix of a com- 
ponent can be obtained for specified port locations. 

When the geometrical shape of a 2-d circuit can be 
considered to be made up of simpler shapes for which the 
Green’s functions are available, the segmentation method [3]»[4] 
can be used for finding the characteristics of the overall 
circuit from those of various components (with simple shapes). 
Details of the analysis using Green's functions and segmenta- 
tion method are discussed in Chapter Two. 

For the analysis of arbitrary shaped 2-d component, nume- 
rical methods are used. Finite element approach [2] and contour 
integral approach [ l] have been discussed in the literature. The 
finite element approach is an extension of finite element methods 
[6] and involves the integration of certain basis functions over 
the entire conducting patch which is divided into numerous sub- 
sections. The contour integral method has been specifically 
proposed for 2-d components and is based on the Green's theorem 
in cylindrical coordinates. The voltage at a point on periphery 
is expressed in terms of line integral along the periphery. The' 
periphery is divided into various sections and the impedance 
matrix is evaluated at these sections. 
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Another approach is the spectral domain analysis [ 7 ], [8] 
which has been used for many planar microwave integrated cir- 
cuits. This method may he used for both the quasi-static and 
fullwave analysis. The quasi-static approach is used to 
compute the capacitance of the structure. The solution of 
the Poisson's equation is obtained in the Pourier transform 
domain. In the full-wave approach, Galerkin's procedure is 
used to deduce determinantal characteristic equation from 
the coupled algebraic equations in the Fourier transform 
domain. 

1.1.2 App lications 

There has been considerable interest in 2-d planar 
components in the last decade and these have found several 
applications in microwave integrated circuits, waveguide 
circuits, ferrite components and microstrip antennas. Some 
of the areas where the 2-d approach has found applications 
are listed below. 

i) : Two-dimensional components can be used for 

the design of high-Q resonators used at microwave frequencies. 
Recently, triangular elements have been proposed for realizing 
resonators and prototype elements for band-pass and band-stop 
filters [ 9 ], [lO]. Gap coupled triangular segments have been 
proposed for filter circuits [ 9 ] and the triangular resonators 
have been used in the design of three-port circulators [ 1 0 ] . 
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ii) Disc ont inuit y Analy sis : Equivalent circuits have been 
reported for most of the discontinuities in strip line and for 
some discontinuities in microstrip line [ll], [ 12] . For 
several discontinuities, equivalent circuits are not available. 
The 2-d approach can be employed to analyze accurately any 
discontinuity. Even when an equivalent circuit is available, 

a more accurate characterization can be obtained by the use of 
2-d approach. It provides a general method for the analysis 
of strip line and microstrip line discontinuities [l3]. 

The concept of cavity type 2-d components, compatible for 
use with waveguide circuits, can be used for the analysis of 
waveguide discontinuities. 

iii) Design__of _Noyel_Gircuit_Configurations : Several new 
designs have emerged from the applications of 2-d circuit 
approach. It has been shown [5] that a circular disc resonator 
can be used as a 3 db hybrid circuit. Also a circular disc 
resonator can behave as a coupled mode filter [l4] when two 
ports are located at right angles and the coupling between two 
modes is obtained by having a very short strip length at an 
open circuited third port. 

iv) *• 2-d circuits, fabri- 

cated on ferrite substrates [ 1 5] , may also be analyzed by the 
2-d approach using the extension of analysis techniques for 
circuits on dielectric substrates. This method has been used 
for obtaining optimum shape of a wideband circulator [ 1 6] . 
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v) Micros trig ^Antennas : Field distribution along the 
periphery of microstrip antennas can he evaluated by treating 
them as 2-d components with magnetic walls. The 2-d approach 
can be used to evaluate the voltages at various points along 
the periphery. The radiation field of these antennas can 
consequently be calculated from the voltage distribution 
[17], [18]. 

1.2 PRESENT INVESTIGATIONS 

The investigations reported in this thesis ewe aimed at 
the extension of Green's function approach for analysis of 
segment shapes not analyzed so far and the use of the 2-d 
analysis in the design and optimization of microwave inte- 
grated circuits. 

Green's functions have hitherto been known for only 
rectangular and circular segments [ 3 ], [ 5 ]. This thesis 
reports the development of Green's functions for some tri- 
angular segments, annular ring segments and some types of 
sectoral segments. Green's functions for the analysis of 
triangular segments are developed in Chapter Three. Three 
types of triangles (30° - 60° right-angled, equilateral, and 
right-angled isosceles) are considered by placing additional 
image sources outside the triangular regions. Green’s fun- 
ctions for circular sectors, annular rings and annular sector 
shaped segments are developed in Chapter Four. 
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Segmentation method for the analysis of 2-d circuits 
involves breaking up of a circuit pattern into segments 
for which G-reen’s functions are known. The scattering 
matrix for the overall circuit is computed from the 
S-matrices of individual components [3], [4]. In this 
method, a considerable effort is spent in computing 
S-matrices for each of the components. These matrices are 
then combined to obtain the overall S-matrix. Considerable 
reduction in computational effort can be achieved if the 
Z-matrices of individual components are combined to give 
the overall Z -matrix from which the overall network 
scattering matrix may be determined. Segmentation method 
that combines Z-matrices of individual components is des- 
cribed in Chapter Pive. 

Two-dimensional resonators are used in MIC's. Resonant 
structures of simple geometrical shapes such as rectangular, 
circular disk and annular ring are used extensively in MIC’s 
[ 19] - [2l]. Planar triangular resonators have been studied 
by Helszajn and James [9J. Spectral domain analysis of 
hexagonal and other microstrip resonant structures has been 
reported [8], [22]. The 2-d approach can be used for the 
analysis of these resonators. The voltages around the 
periphery and input impedances for various modes can be 
obtained. Triangular, rhombic and hexagonal resonators have 
been analyzed using this approach in Chapter Six. 
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Two-dimensional approach has been used to analyze 
strip line and microstrip discontinuities [l3]. The approach 
can be used for the compensation of these discontinuities. 
T-junction discontinuity is studied in Chapter Seven. Nor- 
mally, a rectangular geometry is used at the junction. It 
has been shown that triangular geometry can be used to reduce 
the junction parasitics. 

Stripline and microstrip circuits can be analyzed more 
accurately using 2-d approach. A 3 db branch-line hybrid 
has been analyzed using 2-d approach in [3], [4]. This 
analysis can be used in the optimization procedure to obtain 
improved designs of circuits. Improved designs for 3 db in- 
phase power dividers and 3 db branch-line hybrid have been 
reported in Chapter Eight. 

Chapter Nine- summarizes the new findings reported in 
the thesis and suggests several further extensions of this 


work. 



Chapter Two 


ANALYSIS USING GREEN’S FUNCTIONS AND 
SEGMENTATION 


This chapter presents a review of the Green's function 
approach [l] and the segmentation method [ 3 ], [ 4 ] for the 
analysis of two-dimensional planar circuits. 

2.1 BASIC CONCEPTS 

Consider a two-dimensional circuit shown in Fig. 2.1. 

An arbitrary shaped thin conductor is sandwiched midway 

between two ground conductors spaced 2d apart. There are 

several ports along the periphery. Widths of these ports 

are denoted as W. , W . , ... . The rest of the periphery is 

d 

open-circuited. The xy coordinates are selected parallel 
to and the z-axis is chosen perpendicular to the conducting 
plates. Thus, while the dimensions along x and y coordinates 
are comparable to (or greater than) the wavelength, the 
thickness along the z-direction is negligible. Therefore, 
the fields can be assumed to be constant along the z~ 
direction. 

The general Helmholtz equation, valid for a source free 
region, is given as' ,f 

(V 2 + k 2 )E = 0, k 2 = w 2 pe (2.1) 

where p and e denote the permeability and permittivity of the 

dielectric material and ai is the angular frequency. The above 

^Underlined symbols throughout the thesis indicate vectors or"~ 
matrices. 
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relation holds good for describing the field within the 
dielectric. At the center conductor and at the ground planes 
the tangential components of the electric field are zero. A 
magnetic wall is assumed to exist at the periphery of the 2-d 
component. Since the fields do not vary along the z~direction, 
both E x and E are equal to zero within the dielectric region. 
The fringing field at the periphery is taken into account by 
shifting the magnetic wall by a certain distance from the 
physical periphery. This shift is shown in Fig. 2.2. For 
stripline type planar circuits, the shift 6 in the magnetic 
wall is given by [ll] 


6 = 


2d an 2 
% 


( 2 . 2 ) 


This value of 6 is exact for a semi-infinite patch only, but 
it can be used in most practical cases. For microstrip type 
circuits, the value of 6 is found from the expression for the 
fringing capacitance at an open-ended microstrip [12], 


It may be noted that the fields above and below the 
central patch must be oppositely directed. In the discussions 
below, the E and H can be taken to be the fields on one side 
of the central patch, say the lower side. 


Since E and E are zero, the E-field may be written as 

x y 

I = a z E z U,y) (2.3) 

where a z is a unit vector along the z-direction, and E z is a 
function of x and y coordinates only. Substituting (2.3) in 
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(2.1) and putting 9E z /9z equal to zero, one obtains 

(v^ + k 2 )E z = 0 (2.4) 

2 9 2 9 2 

where vt, = + o’ * Using Maxwell's equation, the mag- 

1 dx 9 y 

netic field can be written as 


V x E 


(2.5) 


which, using (2 .5 ), • reduces to 


— J / z s + a S \ 

3oj tx v ay -x ax — y ' 


( 2 . 6 ) 


where a and a. are unit vectors along the x and the y 
x — y 

coordinates respectively. 

Surface current on a conducting sheet can be obtained from 


the boundary condition 


J s = a x (h 1 - h 2 ) (2.7) 

where, | is a unit vector normal to the boundary and and gg 
are magnetic fields on the two sides of a conducting sheet. 

For the central conductor of a stripline type 2-d circuit, 

= -Hg and thus 

~ 9E 9E 

^ <8^ K + TT %r ) amps/m • (2 - 8) 

For microstrip type planar circuits (after the periphery has 
been extended to take care of the fringing fields), there is no 
magnetic field above the upper conducting plane and so the 
factor of 2 in (2.8) will not appear for obtaining the surface 
current on the conducting patch. 
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The expression for J in (2.8) is valid at all points 
on the central patch including the periphery. For points on 
the periphery, J can be rewritten in terms of components which 
are normal and tangential to the periphery as 
O ^ B ?E 

is = a + IJ- ft) amps/m (2.9) 

where s and n are unit vectors tangential and normal to the 
periphery as shown in Fig. 2.1. For points on the periphery 
where there are no coupling ports, the normal component of 
the surface current must be zero, i.e. 

aE 

= 0. (2.10) 


At a coupling port on the periphery, the planar circuit 
can be excited by a strip line (or micro strip) and the current 
flow (at the coupling port) is normal to the periphery. The 
current flowing in at a coupling port is obtained using (2.9) 
as 


/(fr 


aE, 

j 

an 


•) ds 


( 2 . 11 ) 


where ¥ is the width of the coupling port and ds is the incre- 
mental distance along the periphery* The negative sign in 
(2.11) implies that the current i flows inwards whereas n in 
(2.9) points outwards. 


The characterization of the planar components can be 
carried out in terms of an RF voltage v on the central condu- 
ctor. Since aE z / az =0, v is obtained from (2.3) as 



14 


v = -B d . (2.12) 

Equations (2.4), (2. 10) -and (2.11) can now be written as 
follows : 


(Vj + k 2 )v = 0 


(2.13) 


with 


3y 

an 


= 0 


(2.14) 


for points on the periphery where there are no coupling ports . 
The current flowing in at a coupling port is expressed as 


i = 


/ 


IX 

8n 


ds. 


(2.15) 


w 

Solution to (2.13), with (2.14) and (2.15) as the boundary con- 
ditions, leads to the characterization of stripline type of 
planar components. For microstrip type planar circuits, the 
boundary condition (2.15) is modified so that the factor of 

2 on the right-hand side does not appear. Governing equations 
for waveguide type of planar circuits can be obtained, by using 


similar procedure. 

2,2 GREEN" 1 S PTJNC TI ON" APPROACH 


This approach can be employed when the shape of the 2-d 
patch is relatively simple. The Green's function, which gives 
voltage at any point for a unit source current excitation 
elsewhere, is obtained analytically. When the locations of 
the ports are specified, the impedance matrix of the component 
can be easily derived using the Green’s function. 



15 


2.2.1 Sourc e _ Cur r ents 

If a planar component is excited by a current density J z 
in z-direction at any arbitrary point ( x o »y 0 ) inside the peri- 
phery (as shown in Pig. 2.3), the wave equation can be written 
as 

(v^ + k 2 )v = -jrnp d J z . (2.16) 

When the circuit is excited by a stripline, J denotes a 

Vji 

fictitious RF current density injected normally into the circuit. 

It may be noted that the line current density (J = r~ ), 

ui 3 ^ \ui o n 

being injected into the circuit at coupling ports located on 
the periphery, can equivalently be considered as fed normal 
to the circuit (along the z-direction) with the magnetic wall 
condition 9v/8n = 0 imposed all along the periphery. This 
can be explained as follows. When the stripline is connected to 
coupling port, there is some amount of current in the central 
strip which flows into the planar component. One half of this 
current flows in the opposite direction in each of the two 
ground planes. At the boundary where the strip is connected to 
the planar component, it may be considered that the current 
loops close by equivalent currents along the z-direction. 

Figure 2.4 shows the representation of tangential magnetic field 

(shown as ) at the interface between planar circuit 

and the connecting line by an equivalent current density J g 

(shown as ) in z-direction. These equivalent currents 

may be evaluated as given below. At the periphery, the magnetic 
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field can be written from (2.6) and (2.12) as 


H = — - — (£2, s 
~ jtoiJd v 3n - 


3v 

3S 



(2.17) 


This must equal the magnetic field on the other side of the 

interface between the planar segment and the feeding line. If 

the magnetic wall condition is imposed all along the periphery 

of the planar component, the tangential component of magnetic 

field at the periphery is modelled to be zero. In the present 

model, there is a change in magnetic field at the periphery 

(where coupling ports are located) and the equivalent fictitious 

surface current J in the z-direction, obtained from the boun- 

s 

dary condition (2.7), may be written as 

^ = - S K am P s/m - (2 - 18) 

J is in the negative z-direction in the region below the 
conducting patch as shown in Fig. 2.4. The planar component 
may now be considered to be excited by line currents along 
z-direction at the coupling ports and = 0 can be imposed 
all along the periphery. 

2.2.2 (ir.aaals._EuhatiQua_ahd._Imp a&a&au 

The Green’s function G(r|r 0 ) for (2.16) is obtained by 
applying a unit line current source 6(r - i^) flowing along 
the z-direction in the region below the central patch and 
located at r = r Q . The Green’s function G(x|r o ) is a solution 
of 

(v^, + k 2 )G(r|r 0 ) = -jojpd 6(r - r^) 


(2.19) 
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with, the boundary condition at the periphery given by 

j 

v 

I? - 0 • < 2 - 20 > 

Equations (2.19) and (2.20) yield the Green's functions for 
both stripline type and microstrip type 2-d components. The 
voltage at any point on the planar component may now be 
written as 

v(x,y) = If G(x,y | x Q ,y 0 ) J z (^ 0 t7 Q ) dy Q (2.21 ) 

where J z ( x 0 »y 0 ) denotes a source current density (real or 
fictitious) injected normally and 1 denotes the region of the 
planar component enclosed by magnetic walls. 

When the source current is fed only at ports on the 

periphery, the voltage v at the periphery can be written in 

terms of line current J in z-direction (given by (2.18)) as 

s 

v(s) = - / G(s|s Q ) J s (s Q ) ds Q (2.22) 

C 

where s and s Q are the distances measured along the periphery 
and the integral on the right-hand side is over the entire 
periphery. Since the line current J (s rt ) is present only at 
the coupling ports, equation (2;22) may be written as 

v( s) = - II S(s |s ) J (b ) ds (2.23) 

3 W j 

where the summation on the right-hand side is over all the 
coupling ports and W. indicates the width of the jth coupling 

J 

port. From (2.15) and (2.18), the current i. fed in at the 

J 
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jth port can be written in terms of the equivalent line 
current in the z-direction as 


i-j = -2 / J B (s 0 ) as 0 . (2.24) 

W 3 

If the widths of the coupling ports are assumed to be small 

so that the line current density J is distributed uniformly 

s 

over the width of the port, using (2.24) one gets 


.KM 


for jth port 



Substituting (2.25) in (2.23), v(s) is obtained as 


(2.25) 


v(s) = I p-g 2 / G(s|s ) ds . (2.26) 

j W j W ° ° 

J 

Equation (2.26) gives voltage at any point on the periphery. 
To obtain the voltage v^ at the ith coupling port, average 
voltage over the width of the port is taken, i.e. 


v i = W~ / ds 

1 W. 

# X 

= £ 2 W ' W ^ ^ G (s|s ) ds ds. (2.27) 

^ w i w j mi W. oo 

J J J 

Prom (2.27) above, the elements of the impedance matrix of the 
planar component can be written as 


3-1 


2W . W. 
3- 3 


/ / G(s|s 0 ) ds Q ds. 


W. ¥ 


3 


(2.28) 



20 


The impedance matrix of a component is thus determined using 
the Green's function. The scattering matrix may be evaluated 
from the Z-matrix as [23] 

S = [VY 0 ] (Z - Z Q )(Z + Z^)" 1 [VZ Q ] (2.29) 

where Z Q , [fz ] and [ VT q ] are diagonal matrices with diagonal 
elements given by Z Ql , Z q2 * Z oN ; fZ o1 , fZ o2 » VZ^f 

and l/VZ ol » l/fZ Q2 , *•*» l/VZ oN respectively. Z Ql , Z o2 , ...» 
Z QN represent the normalizing impedances at various ports of 
the network. 


2.2.3 Green '_s_ JTunct, i| ons _f or Jjtectj angular _and. _C1 rcula r 
S egm ents 


Green's functions have hitherto been known for only 
rectangular and circular segments. The Green's function for 
the rectangle shown in Pig. 2.5(a) is given as [3] 


G(x,y |x o ,y o ) 


■itaud 

ab 


CO 00 


cos(k x ) cos (k y ) cos(k x) cos(k y) 

1 1 2 “ — — % ( 2 . 30 ) 


GO «CO 


2 2 2 

k + k ~ k 
x y 


where k x = and k y = The Green's function for the 

circle shown in Fig. 2.5(b) is given as [5] 


G ( P y <J> ! P c y <t> 0 ) 


jojETta 


, . ? v aJ n ( ^ p o )j n (k ttDi p ) cos[n(<j>- * Q )] 

n=o m=1 Tt(a -n /k^ ) (l^-r ) a) 


(2.31) 
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“til 

where J (•) represents Bessel's function of n order , a is 
t for n=0 and 2 otherwise, and k nm satisfy 

5T J n (k nm0>l = °- < 2 -5 2 > 

P =a 

The subscript m in k^ m denotes the mth root of (2.32). For 
zeroth order Bessel's function first root of (2.32) is taken 
to be the nonzero root. 

2 . 3 SEGMENTATI ON METHOD [ 3 ] , [ 4 ] 

The Green's function approach discussed above cannot be 
applied to more general shapes whose Green's functions are 
not known. In many cases, however, the circuit pattern can be 
considered to be made up of simpler shapes whose Green’s 
functions are known. In such cases, the analysis may be 
carried out by using the segmentation method [ 3] ,[ 4 ] . Using 
this method, the continuous connections at the boundaries of 
simpler shapes (which constitute the circuit pattern) are 
replaced by the multiple number of ports covering the entire 
interconnections. The characteristics of circuit obtained 
by this analysis converge to the actual circuit characteristics 
when a large number of ports are considered at the inter- 
connections . 

To illustrate the segmentation method, an example of a 
2-d circuit shown in Fig. 2.6(a) is taken. The circuit shape 
can be considered to be made up of two rectangles. The circuit 
characteristics can be obtained by considering two rectangular 
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components joined at discrete number of ports as shown in 
Fig. 2.6(b). 


The circuit characteristics are obtained by combining 
S-matrices of the simpler shapes joined together. The overall 
S-matrix can be obtained from individual S-matrices by the 
methods used for analyzing multiport circuits. In the multi- 
port connection method [24], the S-matrices of the individual 
components are put together as 


b 

~P 


S 

~PP 

2po ‘ 


a 

~P 

_*c_ 

1 

s „ 

L — cp 

. 


a 

L- — c -» 


where a , b and a , b are normalized wave variables at the 
—p — c ~c 

p external and c internal connected ports. The inter- 
connection constraints are put together as 


b = r a (2.34) 

—ic — — c 

where £ is the connection matrix whose elements are all 

null except the I’s in the entries corresponding to pairs of 

adjacent ports. By eliminating b and a from (2.33) and 

o c 

(2.34), one obtains 


b = S a 

~P -*P ~P 

where 

S = S + S (r - S )~^ S 

^p ^c -cc ; -cp 


(2.35) 

( 2 . 36 ) 


is the scattering matrix of the circuit at its p ports 
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The Green's function approach has been extended to some 
of the triangular and sectoral shapes in Chapters Three and 
Pour. The segmentation method discussed herein is modified 
in Chapter Pive resulting in a procedure where various 
operations performed are in terms of impedance matrices. 



Chapter Three 


GREEN’S FUNCTIONS FOR TRIANGULAR SEGMENTS 

This chapter describes the development of Green's 
functions for some of the triangular segments used in two- 
dimensional planar circuits. For a given shape of 2-d com- 
ponent, derivation of the Green's function requires solution 
of (2.19) with the boundary condition (2.20). Product of 
delta functions 6(x-x Q ) 6(y-y 0 ), cantcuw&i In the right-hand 
side of (2.19), represents a line current flowing along z- 
direction and located at (x Q ,y 0 ). Such sources have been 
called line sources in this chapter. 

In circuits, which are symmetrical about a plane, the 
concept of even and odd modes can be used to simplify the 
analysis of the circuit. For the odd mode half-section, the 
plane of symmetry is replaced by an electric wall at which 
the boundary condition 

G = 0 (3.1) 

is satisfied. Thus, there is need to obtain odd mode Green's 
functions which satisfy (3.1) on the electric wall and (2.20) 
on the remaining sides of the planar segment. 

The method of images has been employed to obtain the 
Green's functions derived in this chapter. 
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3.1 METHOD OP IMAGES [25] 

An analytical solution of the differential equation (2.19) 
is obtained by making the right-hand side a periodic function. 

For this purpose, additional line sources of the type 6(x-x g ) 
•6(y~y_)are placed at points (x_,y_) outside the region of the 
planar component. The additional sources can be thought of as 
obtained by taking multiple images of the line source at (x o ,y Q ) 
with respect to various walls of the planar component. The sign 
of an image source depends on the type of wall across which the 
image is being taken. The sign of the image source is "same for 
an image across a magnetic wall and is reversed for an image 
across an electric wall. The source term in (2.19) gets modified 
and the boundary conditions are satisfied by the voltage v 
produced by the source and its images. The additional sources 
should be all outside the region of the planar component so 
that the solution for G still represents the Green's function 
for the geometrical shape of the 2~d component. The method 
cannot be applied to cases where any of the multiple images 
fall within the region of the planar component itself. 

With these additional images, the source pattern used in 
(2.19) becomes periodic and can be expanded in Fourier series. 

The Green's function can then be expressed as an infinite series 
summation of the functions obtained in the Fourier series expan- 
sion. These are the eigenfunctions of equation (2.13) for the 
given boundary conditions. The coefficients in the series 
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summation for Green's function can then be obtained by- 
substituting these series expansions in (2.19). 

The Green’s functions for some triangular shapes are now- 
obtained using the method of images. 

3.2 GREEN’S FUNCTION FOR 30°-60° RIGHT-ANGLED TRIANGLE 

For the triangular geometry shown in Fig. 3.1(a), the 
Green's function can be obtained by solving (2.19) with the 
boundary conditions given by (2.20). The additional image 
sources are obtained by taking multiple images of the line 
source at (x Q ,y o ) with respect to the three sides of the 
triangle which are magnetic walls. The locations of these 
multiple images are shown in Fig. 3.1 (b). 

It can be seen that the periodicity of the pattern is 3a 
along the y-direction and /3a along the x-direction. Hence, 
ABCDEFGH is chosen as a basic cell which contains 24 line 
sources (all positive*) in 24 triangular regions as shown in 
Fig. 3.1(b). It may be noted that a solution for G with this 
set of multiple sources will satisfy (2.19) in the original 
region of interest. 

For the line source (x 0 ,y Q ) at P in Fig. 3.1(b), the 
coordinates of the image line sources in ODEF can be expressed 
in terms of (x Q »y o ) as 

Since the source at P is positive. 
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for image I 


for image II 


for image III 


for image IV 


for image V 


^ 2a 

= - y 0 + 2 

Coordinates for the other image sources can he directly- 
obtained from these by using symmetry about the axes. 


The Fourier series expressions for each of the 24 line 
sources, which repeat periodically in the two-dimensional 
space, are obtained. The corresponding terms of these 24 
expressions are added together and the resulting expression, 
identical to the right-hand side of (2.19) in the original 
triangular region, can be expressed as 
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8.imid £ 


oi, I eos(k, mx) cos[kp(m+2n)y] 

3 Y3a ni 1 -”* oo n=”* oo 


* {(--l) n cos[k^(m+n)x Q ] cos[k 2 (m-n)y o ] 

+ cos(k^ m x 0 ) cos[k 2 (nH- 2 n)y o ] + 

(-I) 21 *' 11 c os(k^ n x Q ) cos[k 2 (2m+n)y 0 3 } 
where 


k 1 = 


and 


2 it 

f : 5 a 

t - 2 tc 
c 2 3a * 


(3.2) 


(3.3a) 


(3.3b) 


Now substituting *• = -(m+n) , expression (3.2) reduces to 


n . CO 00 

- I I (-1) 

3V3a 


m 


CD — CO 


cos(k^mx) cosEk^Cn-A )y] Tj(x o >y 0 ) 


(3.4) 


where T^(x,y) is defined as 


^(^y) = (~l) £ cosCk^x) cos[k 2 (m-n)y] + 

(**l) m cos(k 1 mx) cos[k 2 (n-Ji )y] + 

(-l) n cos(k^nx) cos[k 2 (S'-m)y] (3.5) 

with the condition that the integers i , m, and n satisfy 


5- + m + n = 0 


(3.6) 


and k^,k 2 are given by (3.3). It can be seen that 
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I I (~1) & cos(k^ Jtx) cosCi^Cm-nJy] ^(x 0 ry 0 ) 

= ll (-l) m cos(k^mx) cos[k 2 (n-£)y] Tj(x 0 ,y 0 ) 

= 1 l (-1 ) n cosCk^nx) cost kg (a -m)y 3 T 1 (x Q ,y o ). (3.7) 

Using (3.7) , expression (3.4) can be rewritten as 


- V*- y) - 


(3.8) 


It can be verified that the functions T^(x,y) satisfy the 
boundary condition, given by (2.20), for the triangle shown 
in Fig. 3.1(a). The functions T^(x,y) are the eigenfunctions 
for this triangle. These eigenfunctions agree with the 
potential functions for an equilateral triangle given by 
Schelkunoff [26]. The Green's function G can be written in 
terms of these eigenfunctions as 


<Mx,y|x ,y ) = l l A ^(x^y) . 

u v rrt r-%— oo v\ co 1 


(3.9) 


m=-“ n=- » 

Substituting (3*9) in the left-hand side of (2.19) and using 
(3.6), one gets 

CO 00 2 

(v2 +k 2 )G = l l [k 2 - (m 2 +mn+n 2 )] A m T. (x,y) . (3. 1 0) 

Since (3.10) and (3.8) are equal for all values of x and y, 
the values of A are obtained by comparison as 


A 

mn 


8j(»pd T^(x o ,y Q ) 


16V”3it 2 (m 2 +mn+n 2 ) - 9V3a 2 k 2 


(3.11) 


Substituting (3.11) in (3.9), one obtains 



G-(x,y |x Q *y 0 ) 


= 8303 pd l l 


T l( X 0* y o ) 

-« 1 Sf^T^ ( m^+mn+n^ ) - 9Y3a^k^ 


(3.12) 


which, is the required Green's function for (2.19) for the 
30°~60° right-angled triangle shown in Pig. 3.1(a). 

3.3 GREEN’S FUNCTION FOR AN EQUILATERAL TRIANGLE 

As in the previous section, the right-hand side of (2.19) 
can he made periodic by placing additional line sources outside 
the equilateral triangle shown in Fig. 5.2(a). The positions 
of additional line sources are obtained by taking multiple 
images of the original line source at (x o »y Q ) with respect to 
magnetic walls at the three sides of the triangle. In this 
case also, the periodicity of the pattern is 3a along the y- 
direction and f3a along the x-direction. The basic cell 
which repeats itself is ABCEEF which contains 12 line sources 
(all positive) in the 12 triangular regions as shown in 
Fig. 3.2(b). 

For the line, source (x Q ,y 0 ) at P in Fig. 3.2(b), the 
coordinates of the image line sources in BCDE can be expressed 
in terms of (x rt ,y^) as 

for image I 


for image II 
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"III 


III 


ia . 12 v , ail 

2 + 2 y o + 4 


3C1 x _ la _ 2a. 

2 o 2 4 


for image III 


"IV 


IV 


- 12. _ n y + S3GL 
2 2 y o + 4 


3d x _ lo 2a 

2 x o 2 4 


for image IV 


"V 




_ x + S3d 

x o + 2 


y o + 2 


2& 


i 


for image V 


Coordinates for the other image sources in the basic cell can 
be directly obtained from these by using symmetry about the 
y-axis. 


Adding the corresponding terms in Fourier series expre- 
ssions for each of the line sources in the basic cell, the 
resulting expression equivalent to the right-hand side of 
(2.19) can be expressed as 


4jj? M 

3V3a 2 


I 

H)—— oo 


£ cos(k^mx) cos [k 2 (m*-2n)y] 


n=- oo 


*{(-1) n cosCk^ (m+n)x 0 3 cos[k 2 (m-n)y 0 3 

+ cos(k,jmx o ) cos[k 2 (mH-2n)y o 3 

+ (-1 ) m+n cos(k 1 nx o ) cos[k 2 (2m+n)y o 3 } + 

cos(k^mx) sin[k 2 (m+2n)y3 { cos(k^mx Q ) sinCls^ (nH-2n)y o 3 
- (-I)™* 11 cos(k 1 nx Q ) sin[k 2 (2m+n)y 0 3 + 

(-1) n cos[k 1 (m+n)x 0 3 sin[k 2 (m-n)y 0 l >] 


(3.13) 
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where k^ and are given by (3.3). Now substituting 
i = -(nn-n) y the above expression can be rewritten as 

_ £2 [(-1 ) m cos(kjinx) cos[k 2 (n-«.)y] T 1 (x Q ,y o ) 

j V jQ* 

+ (-1 ) m cos^mx) sinCk^n-JOy] T 2^o ,y o^J (3.14) 

where T 2 (x,y) is defined as 

^(^y) = (-1 ) 1 cos(k^£,x) sin[k 2 (m-n)y] 

+ (-1) m cos(k^mx) sin[k 2 (n-Ji)y] 

+ (-1 ) n cos (k^nx) sinCk^ £-m)y] . (3.15) 

Integers £ ,m and n in (3.15) may be chosen as in the case of 
T 1 (x 0 ,y 0 ), to satisfy (3.6). It can be seen that 

II(-I) 5 ' cos (k^ £x) Sin[k2 (m->n)y] T 2 (x Q ,y o ) 

= 11 (-1 ) m cos^mx) sin[ (n-£)y] T 2 (x Q ,y o ) 

= 1 1 (-1 ) n cos(k 1 nx) sinCls^ (s. -m)y] T 2 (x Q ,y 0 ) . (3.16) 

Using (3.7) and (3.16), the expression (3.14) reduces to 

- I I[l,(v y 0 > T i (x ’ y) + I 2 (x o’ y o ) T 2 ( x >y)]-(3.17) 

It can be verified that both the functions T^(x,y) and T 2 (x,y) 
satisfy the boundary condition (2.20) for the equilateral 
triangle geometry shown in Pig. 3.2(a). It is seen that 
T.j(x,y) has even symmetry about the x-axis and T 2 (x,y) has odd 
symmetry about the x-axis. Since the line current excitation 
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in the right-hand side of (2.19) is asymmetric, the Green's 
function should contain terms of both T^(x,y) and T 2 (x,y). 


G(x,yJx Qf y 0 ) = l } [A^ T t (x,y) + T 2 (x,y)3. (3.18) 

xa=-^ oo n=~ <» 

Substituting (3.18) in the left-hand side of (2.19), one gets 


(v 2 + k 2 )G = II [k 2 - (m 2 +mn+n 2 )] 

1 9a^ 

* ^ A mn + T^(x,y)] . 


mn 2 


(3.19) 


Since (3.17) and (3.19) are equal for all values of x and y, 

Bmn 

the values of A^ ^are obtained by comparison as 


4j(opd T 1 (x (? ,y 0 ) 


2/2 


16 f 37 t (m +mn+n ) - 9 1 / 3 a k 


(3.20a) 


4iu>nd T 2 (x 0 ,y 0 ) 


1211 16V3ti 2 (m 2 +mn+n' i ) - 9V3a^k‘ :: 

Substituting (3.20) in (3.18), one obtains 


(3.20b) 


G(x,y|x Q ,y 0 ) = 4j«Hd 

“ * T i ( x ,y ) T.(x,y) + T 9 (x ,y ) T 2 (x,y) 

I I Mr 1 — — (3.21 ) 

- ® - -» 16V3tc (m^+mn+n ) - 9V3aK i 

which is the required Green's function for the equilateral 
triangle shown in Fig. 3.2(a). 

Green's Functions m for Even and Odd Mo d e Hal f -Sec tio ns 

First term in the right-hand side of (3.21 ) corresponds 
to the Green's function for 30°-60° right-angled triangle 
given by (3*12). The second term in (3.21) can be identified 
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as one half of the Green’s function for the odd mode half- 
section of the equilateral triangle (i.e. a 30°-60° right' 
angled triangle shown in Pig. 3.1(a) with the side OB being 
an electric wall). This Green's function may be explicitly- 
expressed as 


G(x,yjx o ,y o ) = 8jaipd 


l l 

—00 — 00 


? 2 (x,y o ) T 2 (x,y) 

16V3tc (m +mn+n 2 ) - 9V3a 2 k^ 


( 3 . 22 ) 


3.4 GREEN'S FUNCTIONS FOR AN ISOSCELES RIGHT-ANGLED TRIANGLE 


As in earlier cases, additional line sources are placed 
outside the isosceles right-angled triangle shown in 
Pig. 3.3(a). The positions of additional images are shown 
in Fig. 3.3(b). The basic cell ABCD contains 8 line sources 
(all positive) in 8 triangular regions. Por the line source 
(x ,y ) at P in Pig. 3.3(b), the coordinates of the image 
line source I can be expressed in terms of (x o »y Q ) as 



Coordinates of the other line sources in ABCD can be directly 
obtained from these by using symmetry about the axes. 

Adding the corresponding terms in the Fourier series 
expressions for each of the line sources in the basic cell, 
the resulting expression can be expressed as 




. - 03 

■la»Ud r 
? L 


v „„„ gflrx rniv 

> cos ~ — cos — ■ x - 

u a 


3, 21=— CO xxrr— oo 

where the eigenfunction T(x,y) is given by 


a T < x o 'V 


(3.23) 


m / \ mux mtv , , , o 

T(x,y) = cos p— cos -p- + (-1) 


.nH-n mxx _ nmv 
I cos -p- cos pp-. (3 .24) 


The function T(x,y) satisfies the boundary condition (2.20) 
for the triangle of Fig. 3.3(a). 

As in the previous section, expression (3.23) can be 
simplified to 


_ 22jj& | jf T(x ,y ) T(x,y) . 

2a - ® - co 

As before, 0 can be expressed as 


(3.25) 


CO 00 


G = I I T(x,y) . 


(3.26) 


00 — CO 


The left-hand side of (2.19), then reduces to 


l A V ^ - <?>' - (*»H « x - y ) • 


(3.27) 


The values of A^ are obtained by equating (3.25) and (3.27). 
The Green's function is given by 


. v T(x ,y ) T(x,y) 

G(x,y1x ,y ) = H g o — r 

° o 2[(m 2 +n )it - a 2 k 2 ] 


(3.28) 


which is the required Green's function for an isosceles right- 
angled triangle shown in Fig. 3.3(a). 


A. Green's Function for Odd Mode Half-Section 


The odd mode triangle for the present case has the same 
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shape as shown, in Fig. 3.3(a), with the side OB as an electric 
wall where (3.1) should he satisfied and other sides are mag- 
netic walls where (2.20) should be satisfied. Additional line 
sources are placed by taking multiple images of the original 
line source at (x Q ,y o ) with respect to various walls of the 
triangle. In this case, the basic cell which repeats itself 
is ABCD which contains 32 line sources as shown in Fig. 3.4. 

The signs of the line sources are also shown. Following the 
steps similar to those in the earlier cases, the Green's fun- 
ction for this case is given by 


oo 00 


G(*,y| x 0 ' y o^ = I l 


U(x o ,y Q ) U(x,y) 


u r t 2 2 . 2 . £-1 

oo - oo L (m +n )it -4a k J 
odd(m,n) 


2T 


(3.29) 


where the summation is carried out for odd values of m and n. 


and 


U(x,y) 


COS 


m%x 

2a 


sin 


2a 


hh-el 

-(-1 ) 2 cos 


mix 

2a 


sin 


2a 


(3.30) 


Even mode half -section of an isosceles right-angled triangle 
has the same shape as the original segment. Green's function 
given in (3.28) may be used for this case also. 

B . Green^.s m Function. when^H^gotenuse is jan Jjjlec t ri e _ Wall 

An isosceles right-angled triangle, with the hypotenuse 
as an electric wall and the other two sides being magnetic 
walls, would be needed for odd mode analysis of a square having 
symmetry about one of its diagonals. For example, for the 
square planar segments shown in Fig. 3.5, the odd mode triangle 





octet - i^S| 

; t vimn0% : 
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would be a right-angled isosceles shown in Fig. 3.3(a). For 
this case, the boundary condition (3.1) is valid on the hypote- 
nuse AB and on the other two sides (2.20) is satisfied. The 
positions of the additional sources (alongwith the sign) are 
shown in Fig. 3.6. The basic cell which repeats itself is 
ABCD. By carrying out steps similar to those in the earlier 
cases, the Green's function is found to be 


G(x,y |x 0 ,y 0 ) 


= l l 


Oojixd W(x Q ,y o ) W(x,y) 

2. 2s 2 2, 2 


m=-® n=-» 2L (m +n )it -a k J 


(3.31) 


where 


s mux „ nuy , . sm+n nux „ muy ^ v 

W(x,y) = cos ~ — cos — f — (-1) cos —r~ cos ~~r~ • (3.32) 

d d d d 

3-5 GREEN'S FUNCTION FOR RECTANGLE (ODD MODE) 

A rectangle, with one side as an electric wall and the 
other three sides being magnetic walls, would be needed for 
odd mode analysis of a rectangle having symmetry about an axis 
parallel to and midway between two of its sides. For example, 
for the rectangular planar segment shown in Fig. 3.7(a), the 
odd mode half-section would be a rectangle as shown in 
Fig. 3.7(b). For this case, the boundary condition (3.1) is 
valid on the side 0A and on the other three sides (2.20) is 
satisfied. By carrying out steps similar to those in the 
earlier cases, the Green's function is found to be 




Fig. 3.6 Image sources for isosceles right-angied 
triangle with hypotenuse as electric wall. 
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The method of images used in this chapter is restricted 
to the shapes enclosed by boundaries which are straight lines. 
This is because the only mirror surface which gives a point 
image for a point source is a plane mirror. Only for a few 
polygonal shapes, the images can be uniquely specified with- 
out any -of the image sources falling within the region of the 
planar segment. If all sides of the polygon are magnetic 
walls or all sides are electric walls, the internal angle at 
each vertex should be a submultiple of 180°. For polygons 
with mixed boundary conditions (i.e. some sides as electric 
walls and rest magnetic walls) the angles at all vertices 
where one side is electric wall and the other magnetic wall 
should be 90° or a submultiple of 90°. Other vertices should 
have angles which are submultiples of 180°. Thus, the 
method of images is restricted to rectangular and the triangu- 
lar shapes discussed. 
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In reducing the right-hand side of (2.19) to a periodic 
function, one of the factors obtained, viz. T^(x,y) in (5.8), 
T 2 (x,y) in (3.17), T(x,y) in (3*25), U(x,y) in (3.29), W(x,y) 
in (3.31) and cos(k x x) sin(ky y) in (3*33), is the eigen- 
function which satisfies the boundary conditions. Thus, this 
part of the procedure may be used in evaluation of eigen- 
functions in similar cases. 

The procedure employing the method of images can also be ■ 
used for finding Green's functions for planar segments with all 
short circuit boundaries reported in [2?1 . A similar procedure 
can be used for finding Green’s function for Poisson's equation 
also. 

The Green's functions for some of the triangular planar 
segments developed in this chapter may be used to analyze 
multiport circuits using triangular segments. Other two- 
dimensional shapes which can be segmented into these types of 
triangles (and circular and rectangular geometries) can also 
be analyzed using the segmentation method [ 3 ], [4]. 

These Green's functions have been used for the analysis of 
triangular, rhombic and hexagonal resonators in Chapter Six, and 
for the analysis of triangular T-junctions in Chapters Seven 
and Bight. 



Chapter Pour 


GREEK'S JUNCTIONS FOR CIRCULAR SECTORS, 

AMU LAR RINGS AND ANNULAR SECTORS 

In this chapter. Green’s functions are obtained for 
circular sector, annular ring and annular sector shaped 
segments. Annular ring structures have been proposed for 
use in microstrip antennas [17] and are used in circuits 
like ra trace hybrids. The circular and annular sectors 
are used in microstrip bends as shown in Pig. 4.1, and 
their characterizations are needed to analyze these bends 
accurately. Also, the characterization of an annular sector 
can be used to analyze accurately a tapered line segment as 
shown in Pig. 4.2. The Green's functions of these geometries 
have not been available so far. 

The Green's functions developed in this chapter have 
been obtained by expansion in series of eigenfunctions . 

These Green's functions are valid both for triplate stripline 
type and for open microstrip type circuits. 

4.1 EXPANSION OP GREEN'S FUNCTION IN EIGENFUNCTIONS [25] 

The Green's function G(r|r Q ) is given by the solution of 
(2.19) with the boundary condition (2.20). In this method, 
the Green's function is expanded in terms of known eigen- 
functions of the corresponding Helmholtz's equation given by 
(2.13) and (2.14). 




Fig . 4 ,1 Circular and annular sectors in 
stripitne/riitcrosinp bends - 



F ig . 4 modeled by 
an annular sector 1 
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let represent the normalized eigenfunctions of (2.13), 

2 

which satisfy (2.14), and let k n "be the corresponding eigen- 
values so that 


7 


2 

T 


ip 

n 


+ k! 


n 


ip 

n 


= 0 


(4.1) 


where n represents the two indices defining a particular ip n . 
The normalised eigenfunctions form an orthonormal set, 

such that 

f 1 , if n = m 

// 4 * ’{'m dx dy = < (4.2) 

D 1 0 , otherwise 


where the superscript * denotes complex conjugate, and the 
region Of integration D is bounded by the periphery of the 
planar component, at which ^ n ' s satisfy the boundary condi- 
tion ay / 3n = 0. The Green* s function G satisfies similar 
boundary conditions given by (2.20). Assuming that the 
functions \p form a complete set of orthonormal functions, 
it is possible to expand G(r|r Q ) in a series of let 

G(eIeJ = l A * (r) . (4.3) 

m 

Substituting (4.3) into (2.19), and using (4.1), one gets 

I A (k 2 “ *4 ) ^(r) = Mr - r_) . (4.4) 

m 

Multiplying both sides of the above equation by ^ n (r) and 
integrating over the region D, one obtains 

I A m (k 2 -k^) // if> m (r) ij/*(r) dx dy = -jc^id t},*^) (4.5) 

m D 



49 


which, by virtue of the orthonormal property (4.2), reduces to 


V* 2 - ■ 


* , 


^ n (^) 


(4.6) 


Thus 


A 


n 


3 w Md (Tq) 

Z ? 32 

k n ~ k 


so that 

<M£l£ 0 ) = Dupd 


j ♦ n ( g ) *L ( £q > 


n 


2 2 
4 - r 


(4.7) 


(4.8) 


is the required Green's function esqoansion. For a lossless 
circuit, if> n are real and the complex conjugate is not needed 
in (4.8). 


This method is restricted to cases where the eigenfunctions 
are known. As discussed in the previous chapter, for shapes 
enclosed by straight edges, eigenfunctions can be obtained, 
only if, the internal angle at each vertex is a submultiple 
of %. Eigenfunctions can also be obtained for shapes like 
circles, rings, ellipses, etc. [ 1 73 • 

4.2 CIRCULAR SECTORS 

For the circular sector shown in Fig. 4.3, the set of 
functions which satisfy the boundary conditions is given as 

f vm (p,<t> ) = J v (k vm p )cos v« (4.9a) 

v = rnt/ a (4.9b) 

where J v (*) represents Bessel's function of v order, and 
k vm are chosen to satisfy 
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J v (t vm a) - 0 . (4.10) 

+*Vi 

In above, the subscript m denotes the nr root of (4.10), 

<y ( * ) denotes derivative with respect to the argument, and n 

could be any nonnegative integer. In addition, m can take 

the value zero for n= 0. This corresponds to k =0 and 

v m 

implies that the function has a value of unity throughout the 
sector. The set of functions given by (4.9) are, in general, 
not orthogonal to each other. The functions are orthogonal 
if, and only if, v is an integer which implies that it is an 
integral multiple of a. 

let us consider the case for which a = %/z , where a is 
a positive integer. The set of functions discussed above 
are mutually orthogonal since v takes only integral values. 
To obtain the Green's function using (4.8), these eigen- 
functions must be normalized over the region of the sector. 
The following integral is used for normalizing the eigen- 
function : 

2 

a f a/2,ifm=n^=0 

If J a.( k n.m p > ap = - , - . 

011 \ else i[a 2 -n 2 /k 2 ]j 2 (k a) 

1 n^m n^m 

where n^ = n& , and denotes that v takes only integral 
values. Also 


a 

2 

/ cos n^<j) d<f» 
o 


a , if 
else a/2 . 


= 0 


(4.12) 


I., I < , AlVPljIt 
1. 

v» , 04 

£% . ^ 
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Using (4*1 1) and (4.12), the Green’s function can be written 
as 


g(eIEo) = 


2 a 

K 


V J n. (k n. m p o> J n .< 1 W ) 003 n i+o 008 n i* 

1 X 1 1 X 


to 00 

2 3 ® Ud Z Z p p p p p p 

n=o m=1 %{z - V^.m^Vm ~ k ) ( k ~ ™ a ) 


n . n . m 
x 1 


(4.13) 


The 


angle of the. sector a equals %/s, , and n^ = ns,. The 


eigenvalues ^ m are given by 

V ( v m a > = 0 


and 


n. 


f 1 - 

u. 


if n ± = 0 


if n^ > 0 


(4.14) 


(4.15) 


For A= 1, when the sector becomes a semicircle, the Green's 
function is same as that for the even mode half-section of a 
circle. Thus the Green’s function given by (4.13) is con- 
sistent with the Green’s function for a circle given in (2.31) 


4.3 ANNULAR RINGS 


For the annular ring structure shown in Fig. 4.4, the 
set of mutually orthogonal eigenfunctions, which satisfy the 
boundary condition . (2 . 14) on the periphery, is given by 


pr ^ 

[N^O^a) 


J ' ( k a ) 
n’ nm 




f cos 
i sin 


n<f> 

n4> 


(4.16) 
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where N (■) denotes Neumann's function of order' n, and N'(*) 
n r n 

denotes derivative with respect to the argument. The eigen- 

2 

values are obtained from 




(4.17) 


As in (4.10), n takes nonnegative integral values and m takes 
positive integral values in the above equation also. In 
addition, m can take the value zero for n = 0, and this corres- 
ponds to a unity value of f n throughout the region and 

k =0. 
nm 

To obtain the Green's function using (4.8), the eigen- 
functions described by (4.16) are normalized using 


/ P [^(^ m a) J n (k nffi P) - J^a) \(k m P )] 2 dp = 


£ {(b 2 ~ S*) [N’(k, m a) J (k b) - J'(k a) N (k ,b)]‘ 
* k 2 n ^ 1 ® n nm n' nm n nm 

xim 


- (a 2 - 


n 


-)ClT* (Isu^a) J„(k a) 
' n v nm ' n v nm 


J'(k a) 
n' nm 


TO 


N„(k a)]‘ 
n nm 


> 

(4.18) 


and usual normalizing relations for cos n^ and sin n<J> . The 
Green's function for the annular ring can now be written as 


g (z1e 0 ) 


2 2 

3 a)eu(b -a ) 


+ 


I 1 


cr n (p Q ) ^ ffl (p) cos[n(4>- <J> 0 )3 


n=o m=1 rf(b 2 - n 2 /4)^ m (b) - (a 2 -n 2 /l^)]4(a)3(k 2 m -k 2 ) 

(4.19) 
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where ^^(p) is defined as 


J n< k na p > " J i' k nm a) B n ^nm p > < 4 ' 20 > 

and are chosen to satisfy (4.17). 


To use this Green's function for analyzing an annular 

ring, (4.17) has to he solved repeatedly to obtain k 's. 

nm 

4.4 ANNULAR SECTORS 


The Green's functions for annular sectors (shown in 
Pig. 4.5) can he obtained in the same way as for circular 
sectors. The set of functions which satisfy the boundary 
conditions is given by 


= Cn;(k va a) J (k p) - J' (k um a) K v (k vm P)]oos * 


(4.21) 


where v = nit/a» and k^ m satisfy 


J v <W> JJ < k vm *} 

^ ■ 

These functions are mutually orthogonal if 

a = %/ i , where a is a positive integer. 

can now be written as : 


( 4 . 22 ) 

, and only if, 

The Green's f unc t ion 


o(ElEo) = 

00 CO 

- I I 

n=o m=1 


~2 ' 2 ~ + 2 j Z a) pd 
j toe%( o —ar ) 

CT n ± \m( p 0 ) cos n i *o cOS V 

7t[(b 2 -n^/l^ ia )^ im (b)-(a 2 -n|/^ m )^ m (a)] (^-k 2 ) 

(4.23) 


The angle of the sector a equals %/ s. , and n^ = n£ . The 
function P n m (p) is given by (4.20) and are obtained by 




^ S An annul 
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solving (4.22). For &= 1 , when tho annular sector becomes 
an annular half -ring, the G-reen’s function is same as that for 
even mode half-section of an annular ring. Thus, the Green's 
function for annular sector is consistent with the Green’s 
function for annular ring in the same way as the Green's 
function for a circular sector is with that of a circle. 

4.5 DISCUSSION 

Green's functions for circular sectors, annular rings and 
annular sectors have been obtained by expansion in series of 
known eigenfunctions. Using the Green’s functions given by 
(4.13), (4.19) and (4.23), the Green's function technique 
of analyzing planar circuits and microstrip antennas can be 
extended to shapes incorporating circular sectors, annular 
rings and annular sectors. 

Analysis of circular sector shaped segments using Green's 
functions is only as expensive (computationally) as the ana- 
lysis of circular segments since the eigenvalues in both cases 
are obtained as zeros of the derivatives of Bessel’s functions. 
For annular ring and annular sector shaped segments, the eigen- 
values are obtained by solving the transcedental equations 

(4.17) and (4.22) respectively. Corresponding equation has to 

2 2 

be solved for obtaining each eigenvalue k nm or and thus, 

the procedure is computationally expensive. 



Chapter Five 


SEGMENTATION METHOD USING IMPEDANCE MATRICES 

As discussed in Chapter Two, two-dimensional planar 
components can be analysed by the use of Green’s functions. 
The Green's functions are available for only a few regular 
shapes. The technique of segmentation can be employed to 
analyze planar circuits of shapes which can be considered 
to be made up of those regular shapes. A formulation for 
combining segments to form two-port and four-port circuits 
is given in [ 3 ] and [ 4 ] and could be extended to any general 
n-port network. In this method a considerable effort is 
spent in computing S-ma trices for each of the components. 
These matrices are then combined to obtain the overall S- 
matrix. Considerable reduction in the computational effort 
involved can be achieved if Z-matrices of individual compo- 
nents are combined to obtain the overall Z— matrix from which 
the scattering matrix of the overall network may be deter- 
mined. Segmentation procedures that combine Z-matrices of 
the components are described in this chapter. 

5.1 COMPUTATIONAL EFFORT IN SEGMENTATION USING S-MATRICES 

In the segmentation method using S-matrices, the 
S-matrices of the individual components are obtained and the 
overall S-ma trig is computed using the method discussed in 
Sec. 2.3. The c-ports in equation (2.33) would refer to the 
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ports at the boundaries between the regular shaped components. 

In this method, the solution of (2.36) requires inversion 
of a matrix of order equal to the number of interconnected 
ports. Let us consider an example shown in Fig. 5.1 for illu- 
strating the total computational effort needed in the segmen- 
tation method. This network is a planar circuit version of a 
compensated in-line power divider [28]. Segment B and parts 
and O 2 of the segment C are quarter wave transformers with 
characteristic impedances equal to Z /(2^^) and (2^^)Z Q , 
respectively. Segments A,D and E are portions of outgoing 
transmission lines (with characteristic impedance Z Q ). These 
three segments are considered as planar components in order 
to take into account any higher order modes that may be 
excited by the discontinuities in the line width existing at 
the ends of three transformers. For better accuracy, each 
external port is divided into six subports for obtaining the 
Z-matrix. The six subports are combined together using ideal 
six-way power dividers (not shown in the figure) at each 
external port [ 3 ]. To obtain S-ma trices of individual seg- 
ments, five complex matrices, three of order 12 and one each 
of order 14 and 20, are to be inverted. The number of inter- 
connections in the network is now 44 and so a complex matrix 
of order 88 has to be inverted to obtain the overall scatter- 
ing matrix. 

It has been pointed out by Monaco and Tiberio [293, that 



o 

CM 


X 

o ■ 

CM 

c 





with suitable ordering of rows and columns it is possible . to 
write the connection matrix in the form 



(5.1) 


where I is an identity matrix and 0 a null matrix, each of 

order c/2. They have also pointed out that when two segments 

are being interconnected, submatrix S in (2.33) can be 

~~cc 

written in the block diagonal form as 


See 


M 0 

0 N 


(5.2) 


where M and N are x 2 ma ^ r ices.. In such cases the matrix 
( r ~ S „) ma y be written in the form 


( r - s ) 

- — cc 


-M 


-N 


(5.3) 


The inverse of ( L ~ S ) can now be expressed as 

c c 


( v - s ) 

— -cc 


-1 


1(1 - M N)** 1 (I - £ M)- 1 


,-1 


(I - M N) M(I - £ M) 


-1 


(5.4) 


which requires inverse of two matrices of order ^ and thus 
results in saving of computational effort [29]. 
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This technique can he extended for interconnection of 
more than two segments also, if S can be expressed in the 
block diagonal form of (5.2). This is possible if the net- 
work does not contain any loop with odd number of segments 
in it. The nomenclature ’loop* as used here is a topological 
term meaning any closed path constituted by segments and 
connected subports. Two segments joined together by two 
(or more) subports would constitute a loop with two segments 
in it. More examples of loops are shown in Fig. 5.2. 

In a network, if there is no loop with odd number of 
segments in it then the above technique can be employed to 
reduce computational effort. The circuit shown in Fig. 5.1 
does not contain any loop at all and combination of S- 
ma trices using (5.4) would require inversion of two complex 
matrices each of order 44 only. 

5.2 SEGMENTATION USING Z-MATRICES 

The computational efficiency of segmentation using S~ 
matrices can be improved if the Z-matrices of individual 
components are combined to yield the overall Z-matrix. 

5.2.1 Int ere onnect ion__Using w Z-ma trices 

In a general network of segments, the Z-matrices of 
various components can be written together as 


V 


z z 


i 

~p 


“TP ~PC 


“P 

V 

L ~c J 


_ — cp ^cc_ 


. i c. 


( 5 . 5 ) 
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where v^ and are the voltage and the current variables at 

the p external ports, and v and i are the corresponding 

o c 

variables at the internal connected ports. Equation (5.5) 
does not take into account the constraints imposed by the 
interconnections. The constraints are : (i) the voltages 
at the two connected ports are equal and, (ii) the currents 
are equal and oppositely directed. Eor all the internal 
ports, these constraints can be expressed as 


v c = 0 (5.6a) 

and 

1 2 i c = 0 (5.6b) 

where and matrices with c/2 rows and c columns are 
given by the circuit topology. In these matrices, each row 
describes a connection such that all elements in a row are 
zero except the two corresponding to the two connected ports. 
The two nonzero entries in a row are 1 and -1 for matrix 
. For matrix IV, , both of the nonzero entries are 1 . 

Substituting for v from (5.5) into (5.6), one obtains 


Z i 
— cc — c 


- X, Sop ip' 


(5.7) 


Multiplying (5.6b) by j (= f^-f) and then combining with (5.7), 
the relation between and i can be expressed as 


r z 


r z 

■H — cc 


"i -op 

3 I 2 

0 


(5.8) 



where 0 is a null matrix with c/2 rows and p columns. In 
(5.8), the matrix on the left-hand side is a square matrix 
of order c x c. The current variables i can thus be written 

""C 


as 


1 = 
~c 


-1 


r. z H 

—1 — cc 


i 

i 

Jh-j 

i 

- 3 h _ 


C 




(5.9) 


It should be noted that the Z-matrices for lossless 
circuits are purely imaginary. Equation (5.6b) was multi- 
plied by j so as to make the matrix to be inverted in (5.9) 
purely imaginary. Multiplication and inversion of purely 
imaginary matrices can be carried out with the same computa- 
tional effort as required for real matrices. 

Substituting the expression for i from (5.9) into (5.5) 


Z = Z 
“P ~PP 


Z 

^c 


the overall network impedance matrix is obtained as 

r n Sec i 1 r - r > ^ 

jr 2 o 


(5.10) 


This requires the inversion of a matrix of order equal to the 
number of interconnected ports. For lossless circuits, this 
combination requires inversion of a real matrix of order c x 

5.2.2 Combinatio n U sings Subgrou pe d Interc onn ections 

The computation effort required for combination of Z- 
matrices can be reduced further if the connected ports are 
subgrouped as described below. 
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The 'c' connected ports are divided into groups ’q’ and 


’r' , each containing c/2 ports. This is done in such a way 
that q^ and r^ ports are connected together, q^ and ^ ports 
are connected together and so on. This involves reordering 
of rows and/ or columns for ' .Z „ , Z and Z as given in 
(5-5). The Z-matrices of individual components can now be 
written together as 


v 

— p 


v 

-q 


v_ 


Z ^ Z Z 
— pp — pq — pr 

z z z 

-qp — qq -qr 

Z Z Z 
-rp — rq — rr 


x 

~P 


^q. 


(5.11) 


In this formulation, the interconnection constraints can be 
expressed in a much simpler form as 

= 2 ^ , (5.12a) 

and 


i + i =0 
— q — r — 


(5.12b) 


Substituting (5.12) into (5.11) and eliminating v^, v^ and 


i , one obtains 


or 


^ ^qq ~ ^qr ~ — rq + ^rr^-^q ^^rp ~ ^qp^%> 


iq = ^^qq ~ ^qr ” ^rq + — rr^ ^^Tp ~ ^qp^^p * (5*15) 

Substituting (5.13) into the first equation of (5.11) and 
using i. r = - iq, the Z-matrix of the overall network is 


obtained as 
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%P + ~ %r^%q ■* ^qr ~ ^rq + — rr^ rp ~ ^qp^ * 

(5.14) 

The relation between v^ and i^, obtained by substituting 
(5.13) into the second equation of (5.1 1)» is given as 

^■q ~ ^ ^qp + ^^qq ~ ^qr^^qq ~ ^qr ~ ^rq + ^-rp”^qp * 

(5.15) 

In the above formulation, the order of the matrix to be inver- 
ted is c/2 as compared to the case of (5.10) where a matrix of 
order c needs to be inverted. 


5.2.3 


The two methods of combining Z-matrices are illustrated 
by combining two three-port resistive subnetworks shown in 
Fig. 5-3. The Z-matrices of the two components are 


and 



(5.16a) 


(5.16b) 


Using the first method for interconnection using Z-matrices , 
the component characterizations are put in the form 
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, 

< 

1 


‘4 0 | 3 3 0 0 

, 

’ i 1 ' 

v 5 


! ’ 
0 5 j 0 0 3 3 



1 

<? 1 
ro , 


3 0 1 7 3 0 0 

1 


i 2 

t> 

= 

3 0 j 3 5 0 0 


i 3 

v 4 

! 

o 

o 

o 

U4 


i 4 

- v 6_ 


0 3 J 0 0 3 5 J 


. i 6 J 


Matrices and IV, are obtained as 


L 0 1-1 


-1 


0 


v. 


L v 6 


= 0 


and 


1 0 0 


L o i 1 


■6 1 


= 0 . 


(5.17) 


(5.18) 


(5.19) 


The overall Z-matrix is obtained by substituting (5.17) - 
(5.19) in (5.10). It may be pointed out that the matrix 
need not be multiplied by j since & cQ is real. Thus, 
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‘ 4 

0 “ 


"3 

3 

0 

o' 

z — 








“P 

. 0 

5 _ 


.0 

0 

3 

3_ 


' 7 

3 

-3 

-5 ' 

-1 

' 3 

-3 ~ 

3 

5 

-7 

-3 


3 

-3 

1 

0 

0 

1 


0 

0 

_ 0 

1 

1 

0 _ 


L o 

0_ 


or 


“P 


1 

4 

0 



1 

O 

CD 

— * 

O 

CD 


”3 1 " 

|r| . 



4* 

i 







108 





_0 

5 



108 -108 


1 4 


(5*20) 


This can be verified to be the Z-matrix of the combination. 
For obtaining the Z-matrix of the combination using sub- 
grouped interconnections, the connected ports are divided 
into two groups as discussed earlier. The Z-matrices for 
this purpose can be written together as 
r 


1 


v c 


v. 


L v 3 


4 0 j 3 0 i 0 3 

I I 

0 5)0 3 1 3 0 

» 

3 0 17 0 | 0 3 

0 3 ! 0 7 j 3 0 

0 3 10 3 15 0 

I I 

3 0 j 3 0 10 5 



i 1 


i 5 


*2 


i 4 


± 6 


. i 3„ 


(5.21) 


and the overall Z-matrix is given by substituting various 
submatrices in (5*21) into (5.14),i.e. 




"4 O' 


3 -3 ~ 


"12 -6 “ 

-1 

' -3 3“ 

= 


4- 







- 0 5. 


V*1 

L 


. -6 12 j 


. 3 ~3_ 


’4 0 


-1 1 ' 


’3 r 




= 


_° 5 _ 


_ 1 " 1 


1 4 


( 5 . 22 ) 



71 

-•which, is the Z-matrix of the combination . 

5.2.4 Combination^of^External^Sub^orts 

In many circuits (e.g. in Fig. 5.1), the external ports 
are divided into multiple subports for higher accuracy. The 
combination of Z-ma trices as discussed in Sections 5.2.1 and 
5.2.2 yields the overall Z-matrix at the external subports. 

The multiple subports at each external port are combined as 
follows. It is assumed that only the TEM mode is present in 
the uniform transmission line at the location of the external 
port. This is a valid assumption if the external ports are 
taken at a certain distance away from the planar circuit 
since any higher order modes present at the discontinuity at 
the junction of the transmission line and the planar circuit 
would decay along the uniform transmission line. For con- 
sidering the external ports at a certain distance from the 
discontinuities, the intermediate transmission lines are 
modelled as rectangular planar segments (as shown in Fig. 5.1 
for the power divider circuit ). 

Since only TEM mode is assumed to be present at an exter- 
nal port, the voltages at the subports of an external port are 
same. Thus, one can makes use of the parallel connection of 
subports which implies that the total current injected in a 
port gets divided into its various subports. This combination 
requires inversion of the overall Z-matrix to obtain admittance 
matrix. In general, if ports I and J are divided into subports, 
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I - ...} and J={j^, 32 » } then the term 

T XJ the overall admittance matrix is given as 

Y IJ = l l Y ia (5.23) 

kei s-ej 

where y^ are the terms of the admittance matrix with multi- 
ple subports. The overall scattering matrix may be obtained 
from the admittance matrix Y using [ 23 ] 

2 = [fz 0 ] (E 0 - Ddo + T)' 1 [Vt 0 ] (5.24) 

where [/Z Q ] and [V"Y q 3 have been defined in (2.29) and is 
a diagonal matrix given by Y q1 , Y o2 , ...» Y qN . Y q1 , Y o2 ^ * - Y ohI 
represent the normalizing admittances at various external 
ports. s 

If the external ports are located such that the TEM mode 
approximation is not valid, the voltage and current density 
over the width of an external port can be expanded in terms 
of higher order stripline modes. The higher order stripline 
modes decay away from the junction along the transmission line 
and their characteristic impedances are reactive. Therefore, 
these higher order modes are considered as terminated by these 
reactive impedances at the port location itself [l6]« However, 
for purposes of characterization of the planar circuit, it is 
convenient that the external ports be chosen at the locations 
where only the dominant mode exists, so that the procedure 
described earlier is applicable. 
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5.3 DISCUSSION 

For analyzing the circuit shown in Fig. 5.1 by the inter- 
connection of Z-matrices discussed in Sec. 5.2.1, two real 
matrices of orders52 and 18 and a complex matrix of order 3 
need he inverted. It is seen that there is considerable 
saving in computational effort as compared to segmentation 
using S-matrices. In the method using subgroup ed inter- 
connections discussed in Sec. 5.2.2, inversions of two real 
matrices of orders 26 and 18 and of a complex matrix of order 3 
are required for analyzing the circuit shown in Fig. 5.1. 

Orders of matrices that need to be inverted in different 
methods are compared in Table 5.1. It may be noted that seg- 
mentation using Z-matrices with suitably subgrouped inter- 
connections involves significantly smaller amount of computa- 
tions as compared to the other methods available earlier. 
Moreover, this method can be used for any arbitrary topology 
whereas, as pointed out in Sec. 5.1, the method using S~ 
matrices with subgrouping of connected ports can be used 
only when the network does not contain any loop with odd 
number of segments in it. 

The proposed method using Z-matrices with suitably 
subgrouped interconnections has been used for segmentation in 
Chapters Six to Eight. 
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Table 5.1 Orders of matrices to be inverted in various 
methods of segmentation for the circuit of 
Fig. 5.1 


Method 

Segmentation using 

Segmentation using 


S-ma trices 

Z-ma trices 

Without 

subgrouping 

20 x 20 complex 

52 * 52 real 

the connected 
ports 

3 x ( 12 * 12 ) complex 

18 x 18 real 

14 * 14 complex 

3 x 3 complex 


88 x 88 complex 



With suitable 20 * 20 complex 26 * 26 real 

subgrouping 

of the connected 3 * ( 12 x 12 ) complex 18 * 18 real 

ports 

14 x 14 complex 3 x 3 complex 

2 x ( 44 x 4 . 4 ) complex 


Chapter Six 


TRIANGULAR, RHOMBIC AND HEXAGONAL 
RESONATORS 

Two-dimensional planar resonators are used in 
microwave integrated circuits. These resonators are 
generally of a regular shape. In these resonators, the 
electromagnetic energy is confined in the dielectric 
medium between the ground plane(s) and the planar patch. 
Resonant structures of simple geometrical shapes such as 
rectangular, circular disk and annular ring have been 
extensively used in oscillators, filters and circulators 
[l9] - [2l]. Besides these, resonant structures of other 
geometrical shapes also find applications in microwave 
networks. Availability of different shapes provides 
added flexibility in the design of circuits. Planar tri- 
angular resonators have been studied by Helszajn and 
James [ 9 ] and proposed for use in band-pass and band-stop 
filters. It has been found that a strip line circulator 
using an apex coupled equilateral triangular resonator 
has a bandwidth three times more than that of a circular 
disk circulator [ 1 0 ] . For many applications, the circular 
resonator is not suitable due to anharmonius relationship 
between its mode frequencies [ 30 ]. Hence, there is con- 
siderable interest in the study of resonant structures of 
different shapes [8]. Spectral domain analysis of 
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hexagonal resonators has been reported recently [22]. 

In this chapter, planar resonators of shapes such as 
an equilateral triangle, a 60° - rhombus and a regular 
hexagon are analyzed. Various modes of resonance for 
these resonators have been found from the theoretical 
analysis. Experimental verification of the resonance 
frequencies has been carried out and these have been 
found to be in close agreement with the theoretically 
predicted values. 

6.1 METHOD 0E ANALYSIS 

The resonators discussed in this chapter are analy- 
zed using the Green's function method. The periphery of 
the circuit is divided into various sections (each side 
of an equilateral triangle divided into twelve sections). 
The impedance matrix of a resonator (with various sections 
on the periphery treated as ports) is obtained. For 
rhombic and hexagonal resonators, segmentation method is 
employed to obtain the impedance matrix. The expressions 
for the Green's functions involve series summation from 
-.to to ® . For numerical evaluation of the impedance 
matrix, summation of terms from -50 to +50 is taken and 
is found to be adequate for convergence. Integrals of 
eigenfunctions (appearing in the Green’s function expan- 
sion) over widths of various ports are needed for evalua- 
tion of impedance matrix. These integrals are listed in 
Appendix A. 
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§£S2£§:S9.£-i£_2-d ^Circuits 

The concept of resonance, applicable to 2-d circuits, 
can be evolved from lumped and transmission line resonators. 
In lumped resonators the reactance or the susceptance, 
between two terminals, is zero at the resonance frequency. 

If the resonator is considered lossless,, the impedance is 
zero for series resonance and the admittance is zero for 
shunt resonance. 

A section of uniform transmission line with both ends 

open, resonates when the length of the line is a multiple 

of half wavelength. The input admittance at any arbitrary 

point along the length of the line (Fig. 6,1) can be 

obtained. For the lowest order mode ( X = 2* ), the input 

admittance is zero for all locations of the feed point 

£ 

except when feed point is located in the centre (x = ?r) . 

At this point, the input impedance is zero. This centre 
point corresponds to the zero of the voltage variation 
along the length for the dominant mode. For higher order 
modes also, the input admittance for all locations of the 
feed point is zero except those corresponding to the zeros 
of the corresponding voltage variation. It is seen that 
as the frequency is increased the input admittances at all 
the points (except those corresponding to zeros of the 
voltage variation) cross zero (from negative to positive) 
at resonance. Correspondingly, the input impedances at all 
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these points change sign from positive infinity to 
negative infinity at resonance. 

The concept of resonance described above for trans- 
mission line resonators, can be applied to 2-d resonators 
also. The input impedances, at various sections on the 
periphery, go to infinity at resonance. Just below a 
resonance frequency, the input Reactances at various 
sections have very high positive values and these change 
sign to very high negative values at the resonance 
frequency. This is valid for all sections except when 
mid-point of a section corresponds to zero of the voltage 
distribution. 

The voltage distribution around the periphery at a 
desired frequency of interest can be found from the 
impedance matrix. This requires the location of the feed 
point. If there are more than one feed points, relative 
magnitudes and phases of the current fed have to be known. 

For analyzing resonators of shapes such as a 60°- 
rhombus and a regular hexagon, geometrical symmetry has 
been utilized and the even and odd mode half-sections have 
been analyzed separately. The even mode half-section has 
a magnetic wall at the plane of symmetry and the odd mode 
half-section has an electric wall at the plane of symmetry. 
The impedance matrix is evaluated for the even mode case. 
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For the odd mode case, the impedance matrix is obtained by 
considering the ports along the electric wall to be grounded. 
A general method for finding Z~matrix for a circuit with 
part of the periphery grounded is described below. 

let the ports along the magnetic and electric walls 
be grouped as 1 and 2 respectively. The impedance matrix 
for ports all along the periphery can be partitioned as 


1 

{ 



~ 1 2 


— 

it 

_ —2 

! 

! “21 

-22 


^2 


The subscript 2 denotes ports along the electric wall and 
hence 

y 2 = 0 . (6.2) 

Substituting (6.2) into (6.1), the relation giving the 
impedance matrix for ports along the magnetic walls in the 
odd mode case is obtained as 

3L-j = ~ ~12 ~22 —2 1 ^ if * (6.3) 

The modes of resonance for the three resonator shapes are 
described in subsequent sections. 

6.2 TRIANGULAR RESONATORS 

Planar resonators of equilateral triangle shape have 
been investigated by Helszajn and James [ 9 ]. They have 
studied the field patterns in triangular resonators with 
magnetic wall at the axis of symmetry, i.e. even mode field 
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patterns of an equilateral triangle. For the even mode, 
the voltage along the periphery of the triangle varies 
according to the eigenfunction T^(x,y) given by (3-5). 
The corresponding resonance frequencies are given by [ 9 ] 


k 


m,n,£ 



+ mn + n 


2 


or 

f = — — /m^ + mn + n^ (6.4) 

3afe r 

where c is the velocity of electromagnetic waves in free 
space and 'a' is the side of the triangle. In above, the 
indices £ ,m and n satisfy (3.6). It follows from Chapter 
Three that T 2 (x,y) given by (3.15) is also an eigenfunction 
of the triangle. The voltage variations for these modes 
of resonance are given by (3.15). However, these modes 
are resonant at the same frequencies as the corresponding 
even modes i.e. given by (6.4). It may be noted that 
when m = n, T 2 (x,y) reduces to zero at all points. Thus, 
these odd modes exist only when m ^ n. 


Even modes of resonance are excited when the triangle 
is fed with currents that are symmetrical about the axis 
of sy mm etry as shown in Fig. 6.2(a). Odd modes of resonance 
can be excited by feeding currents which are anti-symmetrical 
about the axis of symmetry as shown in Fig. 6.2(b). These 
odd mode resonances have not been pointed out by earlier 
investigators. At the lowest frequency of resonance, the 
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voltage distributions around the periphery for the even 
and the odd modes of resonance are shown in Fig, 6.3 
(curves 'a* and 1 b* respectively). An arbitrary linear 
scale has been chosen for representing voltages in this 
and in similar figures later. 

For an arbitrary excitation both even and odd modes 
of resonance are excited. The voltage distribution for a 
mixed excitation (-j of the even mode excitation minus 
-jr of the odd mode excitation) is also shown in Fig. 6,3 
(curve 'c')* It is seen that, for this case, the voltage 
variation can be obtained as y of even mode variation 
minus \ of odd mode variation. 

An equilateral triangular resonator has been fabri- 
cated in stripline configuration on polystyrene substrate 
(e r = 2.47, height d = l/l6", side of triangle, A=6.055 cm). 
Photograph shown in Fig. 6.4 includes the lower plate of 
substrate for this resonator. The frequencies of reso- 
nance have been measured experimentally. The resonance 
frequencies are detected using a network analyzer set-up 
as shown in Fig. 6.5. At the resonance frequency high 
fields are excited and because of losses in the resonator, 

, a dip is observed in the reflected signal. To obtain 
theoretical values of resonance frequencies, the magnetic 
wall is considered to be shifted outwards from the physical 
periphery by an amount 6 given by (2.2). The effective 
side of the triangle 'a' is given by 









Triangular, rhombic and hexagonal resonators 
with the top dielectrics and ground planes 
removed . 




2 6 *5 Network analyzer set-up for measuring resonance ' ' | 


frequency . 
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a = A + 2/5 6 ( 6 . 5 ) 

where A is the actual side of the triangle. Theoretical 
and experimental values of resonance frequencies have heen 
compared in Table 6.1. 

6 .3 RHOMBIC RESONATORS 

A 60° — rhombus can be analyzed by segmenting it 
into two equilateral triangles. However, it has been 
analyzed by considering its even and odd mode half-sections. 
The even and odd mode analysis is preferred since the even 
mode half-section corresponds to an equilateral triangle 
whose analysis is already available. Further, this 
approach provides a better insight into the nature of 
modes. 

The even mode half-section of a 60° rhombus is an 
equilateral triangle with all sides as magnetic walls 
whereas the odd mode half-section is an equilateral 
triangle with one side as electric wall and other two 
sides as magnetic walls. These even and odd mode half- 
sections are shown in Fig. 6.6. The rhombus will have all 
the possible modes of resonances of an equilateral trian- 
gle and also those of its odd mode half -section. The 
impedance matrix of the odd mode half— section is obtained 
by using the procedure outlined in Sec. 6.1. 
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Table 6«1 Comparison of theoretically calculated and 
experimentally measured values of resonance 
frequencies for an equilateral triangular 
resonator. 

Effective side a = 6.298 cm, e r = 2.47, d = l/l6 H 


Mode 
m ,n 

A* 

Resonance 

Theoretical 

Frequency 

Exp erimental 

1 ,o 

1.5a 

2.021 GHz 

2.059 GHz 

1 , 1 

af3/2 

3.500 GHz 

3.524 GHz 

2,0 

0.75a 

4.041 GHz 

4.038 GHz 

2,1 

3a/2V5 

5.346 GHz 

5.356 GHz 

3,0 

0.5a 

6.062 GHz 

6.041 GHz 

2,2 

aV3/4 

7.000 GHz 

7.084 GHz 

3,1 

3a/2V3 

7.285 GHz 

7.280 GHz 

4,0 

3a/ 8 

8.082 GHz 

8.047 GHz 


* X = X o /fe r = e/(fVe r ) 
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The modes of resonance corresponding to the even 
mode are same as those of the equilateral triangle. The 
resonance frequencies corresponding to the odd mode are 
obtained separately . The lowest resonance frequency 
corresponds to the odd mode. The voltage distribution 
along the periphery for the lowest order mode is shown in 
Fig. 6.7. The next higher resonance frequency corresponds 
to lowest resonance frequency of the equilateral triangle 
(i.e. even mode half -section) . Since the equilateral 
triangle has two modes (even and odd) at many resonance 
frequencies, the rhombus also has two modes at these 
frequencies. These modes are designated as even-even and 
even-odd, where the first word 'even* refers to the even, 
mode of the rhombus and the second word 'even' or 'odd’ 
refers to the corresponding mode of the triangle. The even- 
odd mode is excited by feeding currents which are anti- 
symmetric about the longer diagonal (xx) and symmetric 
about the smaller diagonal of the rhombus as shown in 
Fig. 6.8(a). The even-even mode can be excited as shown 
in Fig. 6.8(b). The voltage variations around the periphery 
for the even-even and even-odd modes are shown in Fig. 6.9 
(curves 'a' and * b ' respectively). 

The lowest order mode (i.e. odd mode) of the rhombus 
(shown in Fig. 6.7) has quarter cycle variation along each 
side. Hence, this may be called (-£-, i) odd mode. The next 
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higher even-even mode shown in Pig. 6.9 (curve ’a') has 
half cycle variation along each side and can be termed as 
(1,1) even-even mode. The even-odd mode resonant at the 
same frequency has quarter cycle variation along each side 
(Fig. 6.9, curve 1 b ' ) and is termed as (-j, even— odd 

mode. Other higher order modes are also named in the 
same fashion. The positions of zeros around the periphery 
are shown in Fig. 6.10 for various modes. The feed point 
locations have been indicated by arrows. Some higher order 
modes have same number of cycles along each side as another 
lower frequency mode of the same type. For instance, the 
even-even modes at x = 1.5a and X = 0.75a have half 
cycle variation along each side. Though both of these may 
termed as (1,1 ) even-even modes, it is seen that the mode 
at the higher frequency has a non-monotonic voltage varia- 
tion between the positive and the negative peaks along the 
periphery. The variation of voltage along the periphery for 
the even-even mode at *■ = 0.75 a is shown in Fig. 6.11. 

To distinguish it from the lower order (1,1) even-even mode, 
this mode is termed as (1,1)* even-even mode. Similar beha- 
viour has been observed for other higher order modes also, 
where the number of cycles along a side of the rhombus for 
a higher frequency mode is same as that for a similar mode 
at lower frequency because of non-mono tonic variation of 
voltage between positive and negative peaks. In such cases, 
the high frequency mode is termed with a superscript *• 













even™ even mode 
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It is worthwhile here to compare the modes of 
resonanc es of a rhombus with those of a square resonator 
with side a. At the lowest resonant frequency ( x = 2a), 
both (1,0) and (0,1) modes are .resonant in a sq ua re 
resonator. In these modes, the voltage has half cycle 

variation along two parallel sides and is constant along 

* 

the other set of two parallel sides as shown in Fig. 6. 12(a). 

The (1*0) and (0,1) modes can be combined to give the 
(-&•,•£•) modes shown in Fig. 6.12(b). Conversely, each of the 
modes shown in Fig. 6.12(a) can be thought of as obtained 
from the two (i,v) modes in Fig. 6.12(b). In a rhombus, 
the two (-^,-J-) modes, i.e. (i,v) odd mode and (v,y) even- 
odd mode do not resonate at the same frequency and hence 
these cannot be combined to give a (0,1) or a (1,0) mode. 

In a square, the (1,0) and (0,1) modes exist because the 
two (-£-,i!r) modes are resonant at the same frequency (A=2a). In a 
rhombus ,( 1 , 0 ) and (0,1) modes do not exist because the two (i,£) 
modes resonate at frequencies given by x= 2,34a and 
X= 1.5a. 

A 60°-rhombic resonator has been fabricated in strip- 
line configuration on polystyrene substrate (e r = 2.47, 
height d = l/l6 w , side of rhombus, A = 3.755 cm.). The 
photograph shown in Fig. 6.4 includes the lower plate of 
substrate for the rhombic resonator also. The frequencies 
of resonance have been measured experimentally. The effective 
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side of the rhombus is obtained by considering the ' 
magnetic wall to be shifted outwards from the physical 
periphery by an amount 6. Effective side *a* is thus given 
by 


a 


= A + -4 6 

Y3 


( 6 . 6 ) 


where A is the actual side of the rhombus. Theoretical 
and experimental values of the resonance frequencies for 
the rhombic resonator have been compared in Table 6,2. 


6.4 HEXAGONAL RESONATORS 

A regular hexagonal shape has been used as a microstrip 
resonator and as a circulator element C 3 1 3 . Hsu et al.[32] 
have obtained the lowest order resonant frequency of a 
microstrip resonator in terms of effective hexagon side 'a' 

as 

f = 2.011 - • (6.7) 

2-n: a fe r 

Sharma and Hoefer [223 have recently reported analysis of 
regular hexagonal microstrip resonators. They have used 
the quasi-static formulation in the spectral domain to 
determine the capacitance and effective hexagon side. The 
lowest order resonant frequency has then been determined 
using (6.7). 

In this section, the even and odd mode analysis for a 
regular hexagon is presented. The even mode half-section is 
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Table 6.2 Comparison of theoretically calculated and 

experimentally measured resonance frequencies 
for a 60° -rhombic resonator 

Effective side, a = 3.917 cm, e r = 2.47, d = 1/ 1 6" 
Mode X _ Resonanc e freq uency 





Theoretical 

Exp er im en ta 1 

(£>*) ^d 


2.34a 

2.081 

GHz 

2.129 

GHz 

(•£•,£) even-odd 
(1,1) even-even 

} 

1 .5a 

3.249 

GHz 

3.235 

GHz 

(1*,1*) odd 


1 .023a 

4.764 

GHz 

4.765 

GHz 

(1,1) odd 


0.908a 

5.369 

GHz 

5.327 

GHz 

(2,2) even-even 


0.866a 

5.627 

GHz 

5.639 

GHz 

( 1 , 1 )* even-even 
(l-s,1-£-) even-odd 

1 

( 

0.75a 

6.498 

GHz 

6.465 

GHz 

(2*, 2*) odd 


0.66a 

7.383 

GHz 

7.413 

GHz' 

(2,2) odd 


0.596a 

8.172 

GHz 

8.159 

GHz 

(1±,1*)* odd 
(3,3) even-even 
(2^,2-^) even-odd 

) 

0.57a 

0.567a 

8.551 

8.596 

GHz | 

GHz J 

8.549 

GHz 

(2.2) * even-even 

(2.2) even -odd 

I 

0.5a 

9.747 

GHz 

9.737 

GHz 

(3±,3±) odd 


0.482a 

10.109 

GHz 

10.129 

GHz 

(3,3) odd 


0.455a 

10.703 

GHz 

10.704 

GHz 

(4,4) even-even 


0.433a 

11.254 

GHz 

11.194 

GHz 

( , 2 \ )* odd 


0.423a 

11.533 

GHz 

11.526 

GHz 

(3,3)* even-even 


0.416a 

11.714 

GHz 

11.699 

GHz 


(3*, 3*) even-odd 
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a trapezoid which is analyzed by segmenting it in three 
equilateral triangles. The odd mode half-section is also 
a trapezoid with the largest side being an electric wall 
as shown in Pig. 6.13. The odd mode half -section is 
analyzed by finding the impedance matrix for all ports of 
the trapezoid and imposing the electric wall boundary 
condition for the ports along the largest side using the 
procedure outlined in Sec. 6.1. Alternatively, it may be 
analyzed by segmentation method with segment ’a’ (Pig. 6.13) 
being an equilateral triangle and *b' and ’c' equilateral 
triangles with one side as an electric wall. 

At the lowest resonance frequency, both the even and 
the odd modes are resonant. The variations of voltage 
around the periphery for the lowest order modes are shown 
in Pig. 6.14. For many higher order modes also, the odd 
mode resonates at the same frequency as the even mode. The 
positions of zeros around the periphery for various modes 
of resonance studied are shown in Pig. 6.15. For even 
mode resonances, variations along the axis of symmetry are 
also known. These correspond to radial variations (at the 
axis) and are shown in Pig. 6.15 for even modes. The lowest 
order modes of a hexagon are resonant at X = 3.11a. The 
results of Hsu et al. [32] given by (6.7) yield a corres- 
ponding value of 3.124a which is quite close to the value 
obtained by the present method. 
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For circular disk resonators, the even and the odd 
modes resonate at the same frequency* In the case of 
hexagonal resonator , there are certain frequencies where 
either only the even mode or only the odd mode is resonant* 
The nomenclature of modes of the hexagonal resonator is 
based on the modes of a circular resonator* For a circular 
disk resonator of radius ! a f , the modes of resonance are 
given by 


X = 2it/k nm (6.8) 

where k nm satisfy 

J A< k nm a > “ 0 • < 6 ‘9> 

In these modes, m denotes the mth root of (6.9) and 
determines the radial variation of fields. There is n- 
cycle variation around the periphery. 

The lowest order modes of hexagonal resonator have one 
cycle variation around the periphery. For the even mode, 
there is a zero in the center of the hexagon, hence the lowest 
Order modes have been called (1,1) modes. The higher 
order modes have also been named in the same manner. As in 
the case of rhombic resonator, some higher frequency modes 
have similar variations as a low frequency mode. In such 
cases, the high frequency mode is termed with a superscript *. 

It is significant here to compare the modes of resonance 
with those of a circular disk resonator of 


of a hexagon 



107 


radius 'a*. The lowest order (1,1) mode of a circular 
resonator corresponds to x = 3.413a whereas (1,1) mode of 
a hexagon resonates at X = 3.11a. For the (0,1) mode of a 
circle, the odd mode does not exist.- In the present case 
also, only the (0, 1 ) even mode has heen observed. The (4, 1 ) 
mode of a hexagon differs from the corresponding mode of a 
circle so that the voltage in the center of the hexagon is 
nonzero. Resonance frequencies of hexagonal resonator are 
compared with those of a circular disk resonator in 
Table 6.3. 

A regular hexagonal resonator has been fabricated in 
stripline configuration on polystyrene substrate 
(g - 2.47, height d = 1/16", side of hexagon, A = 2.017 cm.). 
The photograph shown in Fig. 6.4 includes the lower plate of 
substrate for the hexagonal resonator also. The frequencies 
of resonance have been measured experimentally. The effective 
side of the hexagon is obtained by considering, as in earlier 
cases, the magnetic wall to be shifted outward by an amount 6, 
Effective side 'a' is given by 

a - A + — 6 (6.10) 

Y3 

where A is the actual side of the hexagon. Theoretical and 
experimental values of the resonance frequencies for the 
hexagonal resonator have been compared in Table 6.4. 
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Table 6.3 


anrcSuLfrlZa^or ? 18 resonanees ° f bezagonal 


Mode 



X for hexagonal 
resonator 

X for circular 
resonator 

(1,1) 

even 

} 



(1,1) 

odd 

3.11a 

3.413a 

(2,1) 

even 

'S 



(2,1) 

odd 

( 

,? 

1.90a 

2.057a 

(0,1) 

even 


1 .5a 

1.640a 

(3,1) 

odd 


1.25a 

1.496a 

(1,2) 

even 


1 .09a 

1.179a 

(3,1)* odd 


1 .09a 

- 

(4,1 ) 

even ^ 




(4,1) 

1 

odd J 

[• 

1 .08a 

1.182a 

(5,1) 

(5,1) 

even *j 

odd J 

i 

0.887a 

0.979a 
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Table 6.4 Comparison of theoretically calculated and 

experimentally measured resonance frequencies 
for a regular hexagonal resonator 

Effective side, a = 2,. 098 cm, e r = 2..47,d=l/l6 ” 


Mode 


Resonan c e F r equ ency 

Theoretical Experimental 


(1,1) 

odd 

j 3.11a 

(1,1) 

even 

(2,1 ) 

even 


l 1.90a 

(2,1) 

odd 


f 

(0,1) 

even 


1.5a 

(3,1) 

odd 


1.25a 

(1,2) 

even 

) 

1.09a 

(3,1)* odd 

J 


(4,1) 

even 

1 

1 .08a 

(4*1) 

odd 

f 


(5,1) 

even 

1 

0.887a 

(5,1) 

odd 

i 


(6,2) 

even 


0.866a 

(3,2) 

even 


0.75a 


2.925 

GHz 

2, 

.903 

GHz 

4.787 

GHz 

4 

.790 

GHz 

6.066 

GHz 

5 

.981 

GHz 

7.272 

GHz 

7 

.203 

GHz 

8.339 

GHz 

8 

.200 

GHz 

8.408 

GHz 

8 

.4-13 

GHz 

10.259 

GHz 

10 

.311 

GHz 

10.506 

GHz 

10 

.400 

GHz 

12.131 

GHz 

12 

.045 

GHz 
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6.5 DISCUSSION 

Lowest order and some higher order inodes of resonances 
of triangular , rhombic and hexagonal resonators have been 
studied. In all the three cases, the measured resonance 
frequencies for various modes are in agreement with the 
theoretically computed values. The difference in the two 
sets of values is mainly due to the tolerance in the 
dielectric constant of the substrate. For ihombic and 
hexagonal resonators, some numerical errors may also be 
present because of finite number of ports considered along 
the interfaces. However, since the number of interconnecting 
ports is twelve along each side of triangular segments, this 
error is expected to be very small. Another source of error 
is that the expression (2.2) used for extension of -the mag- 
netic wall from the physical periphery is not strictly valid 
at the corners. Considering these sources of errors, the 
agreement between theoretical and experimental values of 
resonance frequencies is seen to be good. 

The even mode resonances, of the hexagonal resonator 
considered, correspond to the case of a trapezoidal reso- 
nator. Thus, the results presented could be used for design 
of trapezoidal resonators also. 

For the case of an equilateral triangle, a 60° - rhombus 
and a regular hexagon (the sides being equal in the three 
cases), the value of lowest order resonance frequency is 
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minimum for hexagon and maximum for triangle. For lowest 
order resonance, X = 1.5a for triangle, X = 2.34a for 
rhombus and X = 3.11a for hexagon. 

The voltage variation around the periphery would pass 
through zero an even number of times. For the lowest order 
resonance frequency the voltage would have one cycle varia- 
tion or would pass through the zero twice. The variation of 
A A) even-odd mode for a 60°-rhombus passes through the 
zero twice around the periphery even though it is not the 
lowest order mode. However, this mode cannot be excited 
by a single feed line and two feed points are required to 
excite this mode. The voltage variation for the (0,1) even 
mode of a regular hexagon does not pass through zero around 
the periphery. 

A hexagon has all the resonances of a triangle but only 
the even mode resonances of a rhombus (i.e. same as an 
equilateral triangle). This is because an even number 
of equilateral triangles (i.e. six) form a hexagon whereas 
an odd number of rhombii (i.e. three) make up the hexagonal 
shape. An odd mode of a rhombus cannot be excited in a 
hexagon since it leads to inconsistent boundary conditions 
along AB as shown in Fig. 6.16. For this very reason, all 
the resonances of a triangle are possible in a hexagon also. 

At the lowest resonance frequency of hexagon, both even 
and odd modes are resonant . This behaviour is similar to 
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that of a circular resonator where both even and odd modes 
are resonant (except for (0,m) modes). For some higher 
resonance frequencies, either only the even mode or only 
the odd mode is resonant. This may be because a regular 
hexagon has axial symmetry of 60° only unlike a circle 
which is fully symmetrical. 

The method of analysis used in this chapter may be 
applied to analyze planar resonators of many other shapes 
also. Resonators using annular ring structures, circular 
sectors and annular sectors can also be analyzed by the 
same technique. The Green’s functions obtained in Chapter 
Four may be used for this purpose. 

The resonator configurations studied in this chapter 
could also be used as microstrip antennas. This application 
is, however, not discussed in this thesis. 



Chapter Seven 


T-JUNCTIONS WITH TRIANGULAR GEOMETRY 

T— junctions (or Y— junctions ) are used in stripline and 
microstrip circuits to realize shunt connections of three 
transmission line sections. The T-junctions can be broadly 
considered to be of three different types depending upon the 
impedances of the three transmission lines. The three types 
of T-junctions , as shown in Pig. 7.1, are - i) junction 
where all three impedances are equal, ii) junction with only 
two impedances equal, and iii) junction with all three 
impedances unequal. A rectangular geometry is often used for 
these junctions, like many other discontinuities, the T- 
junction discontinuity also generates parasitic reactances 
which affect the performance of the circuit containing the 
junction. 

Two-dimensional analysis can be used to determine 
accurately the effect of junction reactances and of the 
presence of higher order modes produced thereby. Green's 
function for rectangular geometry given by (2.30) and the 
segmentation method discussed in Sec. 5.2.2 can be used to 
characterize a T-junction. 

It is of interest to know if the use of a triangular 
geometry at the T-junction causes any reduction in the parasitic 
reactances. In this chapter, two kinds of T-junctions with 
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impedance ratios 1 : 1 s 1 and (l/V2):1:1 have been investigated 
for this purpose. For T-junctions with impedance ratio 1 : 1 s 1 , 
it has been found that an equilateral triangular geometry at 
the junction yields the best performance. A right-angled 
isosce3.es triangular geometry yields the best perfor ma nce 
for T-junctions with impedance ratio (l/Y2):1:1. 

2-matrices for the geometrical shapes are obtained 
using the Green's functions given in Chapters Two and Three. 
Integrals of eigenfunctions (appearing in Green's function 
expansions) over widths of various ports are used in the 
evaluation of Z-matrices. These integrals are listed in 
Appendix A. Portions of outgoing transmission lines from 
the junction are also modelled as planar rectangular compo- 
nents. As discussed in Chapter Five^this process accounts for any 
higher order modes excited by the discontinuities. The 
Z-matrices of various components are combined using the 
segmentation method. For higher accuracy, the ports corres- 
ponding to the outgoing line sections are divided into multi- 
ple subports. These subports are combined using the technique 
described in Sec. 5.2.4. The S-matrix is obtained from the 
overall Y -matrix at external ports using (5.24). 

7.1 T-JUNCTIONS WITH EQUAL LINE IMPEDANCES 

In this section, the case of T-junctions with equal line 
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impedances are considered. One of the applications of this 
type of junctions is the circuit for SPDT switches [ 35 ]. 'In 
these junctions, the widths of the three outgoing lines are 
equal and a square geometry is used at the junction. Under 
ideal conditions, the input reflection coefficient at the 
junction should be -1/3 and the transmission coefficient from 
one branch to another should be 2/3. However, due to dis- 
continuity reactances, the performance would differ from 
that of the ideal junction. 

2-d analysis has been carried out for the rectangular 
T-junction. The circuit is divided into segments as shown in 
Pig. 7.2(a). The width of 50 8 line is divided into six 
subports as shown in Pig. 7.2(b). It is found that the refle- 
ction coefficients differ from the ideal value of -1/3 and 
these deviations increase with frequency. The reflection 
coefficient at port 1 is shown as a function of frequency 
in Pig. 7.3 (curve 'a'). 

Since, the widths of the three outgoing lines are same, 
use of equilateral triangular geometry at the junction is 
explored. An equilateral triangle is used at the junction 
and the resulting T- or Y- junction (shown in Fig. 7.2(c)) is 
also analyzed using the 2-d approach. The reflection 
coefficient (equal at all three ports) is very close to the 
desired value of — 1/3 and deviates only slightly at high fre- 
quencies. This reflection coefficient is shown as a function 
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of frequency in Pig. 7.3 (curve •}>•). It is seen that the 
performance of equilateral triangular 1-junction is much 
better than that of a rectangular T-junction. 

To verify the improvement obtained by the use of equi- 
lateral triangle in T-junctions with impedance ratio 1:1:1, 
two T-junctions have been fabricated: one using a rectangular 
T-junction and the other using an equilateral triangular T- 
junction. These circuits have been fabricated in stripline 
configuration on polystyrene substrate (e - 2.47, height 
d = l/l6"). As in the previous chapter, magnetic wall is 
considered to be shifted out from the physical periphery by 
an amount 6 given by (2.2). For 50 lines, the effective 
side of the equilateral triangle is 0.381 cm. and the corres- 
ponding physical dimension is 0.241 cm. In the case of 
rectangular T-junction, these dimensions correspond to the 
square region at the junction. A photograph of the lower 
plates of substrate for the two circuits is shown in Pig. 7. 4. 

For measurement of input VSWR, ports 2 and 3 of T-junction 
are terminated in matched loads and input reflection coefficient 
is measured at port 1 . The variations of VSWR with frequency 
for the two circuits are shown in Fig. 7.5* The theoretically 
predicted variations are also shown in the figure. 

7.2 T-JUNCTIONS WITH IMPEDANCE RATIO (l/f2)1:1 

In this section, T-junctions with one line impedance 






Pig. 7.4 Experimental T-junctions (impedance ratio 1:1:1) 
using rectangular and triangular geometries with 
the top dielectrics and ground planes removed. 




■i « 
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equal to t/lf2 times the other two (equal) line impedances 
are considered* These junctions occur frequently in circuits 
such as power dividers, hybrids, etc. [28]. In these jun- 
ctions the effective widths (in the planar waveguide model) 
of the three lines are in the ratio f2: 1 : 1 . Usually, a 
rectangular geometry is used at the junction. In ideal case, 
the input reflection coefficient at the line with lower 
impedance should be (1 - f2)/(i + ][2) or -0.172. The chara- 
cteristics of a practical T-junction would however be 
different from the desired ideal values because of the dis- 
continuity effects. 

The rectangular junction has been analyzed using 2-d 
approach. The circuit is divided into segments as shown in 
Pig, 7.6(a). As before, for higher accuracy, the ports on 
the 50 n main line are divided into six subports and the 
port on 35.35« line is divided into eight subports. It is 
found that the performance of the junction differs from the 
ideal case and it deteorates with increase in frequency. The 
reflection coefficient at port 1 is shown as a function of 
frequency in Pig. 7.7 (curve 'a'). 

Since the effective widths of the three outgoing lines 
are in the ratio f2:1:1, use of a right-angled isosceles 
triangle (ratio of sides being V2:1tl) is explored. A right- 
angled isosceles triangle (in the effective magnetic wall 
model) is used at the junction. The T-junction (shown in 
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Fig. 7.6(b)) is analyzed using the 2-d approach. The per- 
formance is very close to the ideal behaviour upto very high 
frequencies and the deviation from the ideal case is negli- 
gible. The reflection coefficient at port t is shown as 
a function of frequency in Fig. 7.7 (curve 'b'). It is seen 
that the performance of right-angled isosceles triangular T- 
j unction is much better than that of a rectangular T— junction. 

It should be noted that, for a T— junction with impedance 
ratio (l/V*2):1:1, a right-angled isosceles triangle is 
proposed in the effective dimensions and not in the physical 
dimensions. Only the effective widths of the three outgoing 
lines are in the ratio (l/f 2 ): 1:1 and the actual geometry at 
the junction is not a right-angled triangle. In a typical 
case, the physical periphery and the magnetic wall location 
are shown in Fig. 7. 1 8. 

To verify (experimentally) the improvement caused by 
using a triangular junction, a quarter wave matching section 
is incorporated at the input end as shown in Fig. 7.9. Two 
circuits, one using a rectangular junction and the other 
using a right-angled isosceles triangular junction, have been 
fabricated in strip line configuration on polystyrene substrate 
(e = 2.47, height d = l/l 6 "). As before, the magnetic wall is 
considered to be shifted out from the physical periphery by 
an amount 6 given by (2.2). Various effective and actual 
dimensions of the two circuits are as follows : 













1 29 



Effective 

Actual 

a ) Ree tagjgular^T-june tion 



i) 50^ side of rectangle 

0.381 cm. 

0.241 cm. 

ii) 35.35ft side of rectangle 

0.539 cm. 

0.398 cm. 

iii) Length of rectangular 

segment on the input side 

1 .945 cm. 

1 .945 cm. 

b ) Tr iangular JD; junc t ion 



i) 50 n side of triangle 

0.381 cm. 

0.244 cm. 
(side AS 
Fig. 7-8) 

ii) 35.35ft side of triangle 

0.539 cm. 

0.398 cm. 

iii) Length of the rectangular 

1.910 cm. 

1 .950 cm. 


segment on the input side 

A photograph of the lower plates of substrate for the two cir- 
cuits is shown in Fig. 7.10. 2-d analysis is used to obtain 

the theoretical results. Matching section, geometry used at 
the junction, and portions of outgoing transmission lines are 
modelled as planar segments. The theoretical and experimental 
variations of input VSWR with frequency for the circuit using 
rectangular junction are shown in Fig. 7.11* Experimentally 
observed variation of input VSWR for the circuit using tri- 
angular junction is compared with the theoretical variation in 
Fig . 7.12. 

7.3 APPLICATIONS OF TRIANGULAR JUNCTIONS 

The T— junctions discussed in the earlier sections are 
widely used in various microwave circuits. The performance 




Fig. 7.10 Experimental circuits of (l/V2);1:1 

T-junction with matching section on input 
side using rectangular and right-angled 
isosceles triangular junctions with the top 
dielectrics and ground planes removed. 




heoreticat 
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of such circuits can be improved by using suitable triangular 
shaped junctions. In this section, two examples of these 
circuits are being given to bring out the improvement caused by 
the use of triangular junctions. 

Consider an in-phase power divider circuit using single 
matching sections (Z = Z Q f2) on the output sides as shown in 
Fig. 7. 13(a)* The circuit uses a^ T-junction with impedance 
ratio ( 1/V2 ) : 1 : 1 and hence a right-angled isosceles triangle 
can be used at the junction as shown in Fig. 7.15(h). The 
2-d approach is used to analyze both these circuits designed 
at 6 GHz center frequency. The variations of input VS'WB. with 
frequency are shown in Fig. 7.14 (curve 'a' for the circuit 
using rectangular junction and curve ’b' for circuit using 
triangular junction). It is seen that the center frequency is 
shifted from the design objective in both the cases. Further, 
the circuit using triangular junction has a lower value of 
VSWR (= 1.007) at its center frequency than for the circuit 
using rectangular junction (where VSWR = 1.052). 

Another example of a circuit using T-junctions with, 
impedance ratio (l/V"2):1:1 is a branch-line hybrid shown in 
Fig. 7.15(a). Triangular junctions can be used in this case 
also by making the quarter wave sections curved as shown in 
Fig. 7.15(b). For carrying out its 2-d analysis, the circuit 
should be divided into rectangular, triangular, and annular 
sector shaped segments. However, the characterization of an 










VSWR 
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annular sector shaped segment can be approximated by modelling 
it as a rectangular segment . 2-d analysis for the circuit 

shown in Pig. 7.15(b) has been carried out at 4 GHz. center 
frequency and by modelling the annular sector shaped segments 
also as rectangular segments. The input VSWR’s (equal at 
all ports) for the two circuits shown in Fig. 7.15 are plotted 
as a function of frequency in Fig. 7.16 (curve ’a’ for the 
circuit using rectangular junctions and curve * b * for the 
circuit using triangular junctions). As in the ease of 
power divider circuits, the center frequency of hybrid is 
shifted in both the cases. However, the shift in case of the 
circuit using triangular junctions is much lower (= 220 MHz) 
than that ( = 870 MHz) for the circuit using rectangular 
junctions. The circuit using triangular junctions has a 
lower value of VSWR (= 1.013) at its center frequency than 
for the circuit using rectangular junctions (where VSWR=1.093). 

7.4 DISCUSblON 

It has been shown in this chapter that, the affect of 
parasitic reactances can be reduced for two kinds of T- 
junctions by using triangular geometries at the junction. 

For a T-junction with equal line impedances, an equilateral 
triangle is proposed and for a T-junction with impedance 
ratio ( l/V"2 ) : 1 : 1 , a right-angled isosceles triangle (in 
magnetic wall model) is proposed. Improvement using trian- 
gular shaped junctions has been verified experimentally for 

both the cases. 
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Xi# is perhaps significant to compare the areas of 
junctions with different geometries* for T-junctions with 
equal impedances, effective area at the junction is W 2 for 
the rectangular junction and W 2 /3/4 for the equilateral 
triangular junction where ¥ is the effective width of the 
outgoing transmission lines* Por T— junctions with impedance 
ratio (l/V2):1:l, effective area at the junction is W 2 ^ 
for the rectangular junction and W 2 /2 for the right-angled 
isosceles triangular junction where W is the effective width 
of the two symmetric outgoing transmission lines. Thus, it is 
observed that the reduction in parasitic reactances obtained 
by the use of triangular junctions is associated with the 
reduced junction area. Also, for the two types of junctions 
studied, the corresponding triangles have the least area at 
the junction. To support the argument given above for impro- 
vement in paraeitics, optimization was carried out for T- 
junctions with impedance ratio (l//2):1:1. Three geometries - 
i) rectangular, ii) equilateral triangular, and iii) right- 
angled isosceles triangular were tried at iiie junction (as 
shown in Pig . 7*17)» and their dimensions were treated as 
variable parameters. These three junctions were optimized to 
minimize the deviation of reflection coefficient from its ideal 
value of -0.172. least square optimization procedure described 
in Appendix B is used. In each of the three cases, the minimum 
value of deviation occurs when the dimensions at the junction 
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are minimum for that geometry (i.e. sides of rectangle are 
W 1 and W 2 f side of equilateral triangle is Wj and the longer 
side of the right-angled isosceles tri ang le is. ) . Amongst 

the three optimum values of reflection coefficient, the value 
of input reflection coefficient for the case of equilateral 
triangle lies between the other two optimum values* For a 
typical case (s r = 2*55, d = 1/16", Z 2 = 50 o f = io GHz) 
the input reflection coefficients for rectangular, equilateral 
triangular and right-angled isosceles triangular T-junctions 
have the magnitudes 0.386, 0.214 and 0.189 respectively.- Area 
of the equilateral triangular junction is W^V”3/2 as compared 
to w|/2 and Af 2 for the rectangular and right-angled isos- 
celes triangular junctions respectively . In view of the 
observations summarized above, it may be stated that the equi- 
lateral triangle and the right-angled isosceles triangle 
have the optimal geometry for the T-junctions with equal line 
impedances and line impedances in ratio (l/]f2):1:1 respectively. 
These observations can be extended to T-junctions with arbitrary ;■ 
impedances also. For a T-junction with arbitrary impedances, , 
a triangle with sides equal to the effective widths of the 
three outgoing lines has the least area. It is expected that 
use of this type of triangles would minimize parasitic react- 
ances for T-junctions with arbitrary impedances also. Ana- 
lysis of these junctions can be carried out using contour 

| 

integral method [iJ of 2-d approach. : 

t ‘ 



142 


In both the circuit examples considered in Sec; 7; 3, 
a shift in the center frequency is observed; To obtain 
minimum VSWR at the center frequency, changes in lengths of 
matching sections are needed. Such changes are required 
since in 2-d analysis the matching sections are modelled more 
accurately and the actual dimensions may no longer remain x/4; 
These changes can be obtained by optimizing the lengths of 
various matching sections so that the magnitude of is 
minimized. For the power divider circuit, there would be 
only one variable parameter in optimization; Thus, the 
circuit analyzed would correspond to an optimum design at some 
other center frequency and the VSWR obtained at that frequency 
is the minimum obtainable for that center frequency. However, 
for the branch-line hybrid circuit, there would be two varia- 
ble parameters in optimization; Thus, the VSWR obtained at 
the center frequency need not correspond to the minimum obtai- 
nable for that center frequency. 

Besides power dividers and branch-line hybrids, the T- 
junctions discussed in this chapter are used in stubs, SPDT 
switches, ratrace hybrids and branch-line hybrids. Triangular 
T-junctions can be used in these circuits to reduce the effect 
of parasitics. Improved designs for power divider circuits and 
branch-line hybrid using triangular T-junctions have been 
obtained in Chapter Bight, 



Chapter Bight 

OPTIMIZATION OP STRIPLINE AND MICROSTRIP CIRCUITS 
USING 2-D APPROACH 

Stripline circuits contain many discontinuities whose 
reactances lead to degradation of circuit performance# The 
equivalent circuits for discontinuities are not very accurate 
at high frequencies or when the field perturbations caused by 
the neighbouring discontinuities interact with each other. In 
the previous chapter, 2-d approach has been used to characterise 
T-junctions and for analysis of some circuits which contain 
T-junctions. It has been pointed out there that in the design 
of power dividers and branch-line hybrids the quarter wave 
matching sections do not yield the lowest VSWR at the center 
frequency and optimization needs to be carried out. 

In this chapter, 2-d analysis has been used in the optimi- 
zation of in-phase 3 db power dividers and branch-line hybrid 

circuits [28]. 

8.1 IN-PHASE POWER DIVIDERS 

Three types of power divider circuits have been considered 
as shown in Fig. 8.1. These are : i) a circuit with single 
matching section (Z * Z Q /f2) on the input side, ii) a circuit 
with single matching sections (Z = Z 0 V2) on the output sides, 
and iii) a circuit with a matching section (Z = ^ 0 /2 ) on 

the input side and matching sections (Z = Z Q 2 ) on the 
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output sides. Theoretically (i.e. in absence of discontinuity 
reactances and any higher order modes) each of the matching 
sections should be quarter wavelength long. Output VSWR and 
the isolation between two output ports are improved by 
connecting a resistance of value 2 Z q between the two outgoing 
lines at quarter wavelength away from the junction [28] . 

Bach of the three types of power divider circuits 
(Pig. 8.1) uses a T-junction with impedance ratio ( 1 / ) : 1 ; 1 . 
Thus, their designs can be improved by using right-angled 
isosceles triangular junctions as discussed in the previous 
chapter. Optimization procedure is used to obtain lengths of 
various matching sections so that the input reflection coeffi- 
cient is minimized at the center frequency. Por this purpose, 
matching line sections and portions of line lengths on the 
three outgoing transmission lines are modelled as rectangular 
planar segments. Por obtaining the input reflection coefficien’ 
the symmetry in the circuit is exploited and only the even mode 1 
half -section is analyzed. The case of power divider with mat- 
ching section on the input as well as the output sides, and its 
even mode half-section are shown in Pig. 8.2. Overall input 
reflection coefficient is same as the 1 for the even mode 
ha If -section. A least square optimization method is used to 
minimize |S^| for the even mode with 10 GHz center frequency. 
The optimization procedure used is described in Appendix B. 

The design and performance data (input VSWR's) for the three 
types of circuits is given in Table 8,1. 




F,g 8 2 p ower divider with matching section^ 

at both sides and its 1 

■ ' ( fe) ; ; m 0.d^:|ft ; ci|| | o n £ 
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Table 8.1 Design and performance data of power dividers 
using triangular T-junctions 


(e r = 

2.55, h 

1/16”) 


C onf iguration 

lS 11 J at 
10 GHz 

Bandwidth 
(VSWR < 

1.05) 

Optimized 

lengths of matching 
sections 

Single section 
at input 

0.0043 

9.5 to 

10.5 GHz 

0.425 cm. 

Single sections 
at outputs 

0.0017 

9.55 to 

10.55 GHz 

0.415 cm. 

Matching sections 
at input and at 
outputs 

0.0034 

8.15 to 

11 .3 GHz 

0.436 cm. 

0.429 cm. 
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To reduce output VSWR, an isolation resistance is 
connected between the two outgoing lines from the T-junction. 
This is shown in Pig. 8.3 for the circuit with single matching 
sections at the outputs and its odd mode half-section. The 
lengths of various matching sections (having been determined 
to minimize |S^| ) are now fixed and objective is to minimize 
I S 22 I allowing the isolation resistance and its position 
to be varied. 

It should be noted here that the isolation resistance 
does not affect the performance of the even mode half-section 
since In the even mode case the isolation resistance is open 
circuited in the middle.. Also, to minimize the output VSWR 
it would be sufficient to minimize the t ^^2 ^ f° r the 0< ^ mode 
case since Sgg for the even mode half-section is almost 
zero at the center frequency. 

The analysis of odd mode half-section requires characteri- 
zation of a segment with a resistance connected to ground from 
a specified location. This can be obtained by having a port 
of small width (equivalent to the width of the chip resistor 
to be used) at the position of the isolation resistance and 
terminating this port to ground through a resistance R/2. 

Width of the chip resistor has been assumed to be 0.5 mm in 
the computations. The impedance matrix at the remaining 
ports (with one port connected to ground through a resistance 
R/2) is computed as follows. 




(a) 


Fig 8.3 (a) Power divider with matching sections 

at., outputs;; an^ 
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as 


Let the impedance matrix of the 


segment be partitioned 



subscript r denotes the port to be terminated and subscript 
1 denotes all other ports of the segment. Using 


v _ R . 
r 2 ■ L r 


( 8 . 2 ) 

the impedance matrix at other ports can be obtained as 


= 


* 1.1 


1 


2 + | *1r *ri 

rr 2 


i 


1 


(8.3) 


Various segments are combined using segmentation method and 

optimization procedure is used to minimize |s 22 l for the odd 

mode. The data for the design of isolation resistance for the 

three types of power divider circuits studied is given in 
Table 8.2. 

The frequency variations for the input and output VSWR's 
for the three circuits are shown in Pig. 8.4. 


The experimental verification of the procedure is 
carried out by conducting measurements for a power .divider 
with single matching section on input side. Optimization of 
a circuit (for stripline configuration) is carried out at 
f Q - 2GHz. (d = 1/4”, e r = 10). The concept of scaling 
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Table 8.2 Design and performance with isolation resistance 


Configuration 

Isolation 

resistance 

Position 
(distance 
from jun- 
ction) 

|S 22 | at 
10 GHz 

Minimum value 
of jS 2g | 

(frequency) 

Single 
section at 
input 

80 . 93 0 

0.512 cm. 

0.167 

0.166 (10.1 GHz) 

Single 
section at 
outputs 

134. 97Q 

0.390 cm. 

0.1923 

0.1763(7.2 GHz) 

Matching 
sections at 
input and 
outputs 

131 .36 M 

0.402 cm. 

0.1852 

0.167 (13 GHz) 
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described in Appendix C is used to obtain from this the 
design for microstrip configuration with f = 1 GHz. 

(d = l/4", e r = 10). This circuit design, used for fabrica- 
tion, is a scaled model (Appendix C) of a circuit on 

0,025" thick alumina ( e = 10) substrate with f = 10 GHz. 

r o 

The dimensions obtained from scaling are for magnetic wall 
model (length of 35.35ft segment = 2.194 cm.). The physical 
dimensions of the circuit are obtained using the procedure 
described in Appendix C. The dimensions of the power divider 
circuit (in microstrip conf iguration are shown in Fig. 8.5 
and an outline of the circuit layout is shown in Fig. 8.6. 
The measurements have been carried out using an automatic 
network analyzer*. In these measurements, no compensation 
has been made for connector imperfections. The theoretical 
and measured values of input VSWR are shown as a function of 
frequency in Fig. 8.7. Theoretical and measured values of 
coupling (i.e. J ^ l ) are plotted in Fig. 8.8 as a function 
of frequency. Dotted lines in these figures show maximum and 
minimum values considering the uncertainity in reflection 
coefficient to be 0.06 and in transmission coefficient to be 
0.045* With this uncertainity taken into account, the 
agreement between the theoretically predicted values and 
experimental results is reasonable. 

* Courtesy Mikrowelien-iiab.oratorium, Eidgeridssische I'echnische 
Hochschule ZUrlch. ' ? 
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8.2 BRANCH-LINE HYBRID 

Branch-line hybrids are commonly used in microwave 
circuits. At high frequencies the performance of the circuit 
is limited by discontinuities. This circuit uses T-junctions 
with impedance ratio (l/V"2):1:1. As shown in the previous 
chapter, use of right-angled isosceles triangular geometry at 
these junctions reduces the parasitics. Por using triangular 
junctions in a hybrid, two of the quarter wave sections have 
to be curved as shown in Pig. 7* 15(b). The performance of 
branch-line hybrids using sections of exactly \/ 4 length was 
evaluated in Sec. 7.3. It was observed that the center 
frequency is shifted and the minimum VSWR obtained is about 
1.05. 

The lengths of the two quarter wave sections can be 
modified to obtain optimum performance. All the four quarter 
wave sections are modelled as rectangular planar segments. 
Higher accuracy can be obtained by modelling the two curved 
segments as annular sectors with a 90° sectoral angle. The 
Green's function for annular sectors given in Chapter Pour 
may be used for this purpose. However, modelling of annular 
segments requires equation (4.22) to be solved for each eigen- 
value and thus the procedure is computationally expensive. Por 
this reason, it has not been employed in the work reported here 

Por optimizing the design of the hybrid, there are two 
designable parameters - i) length of the Z q /^2 section, and 
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ii) length of the quarter wave section with impedance Z Q . 

The objectives for optimization are to minimize VSWR and to 
equalize power division at the center frequency. There are 
two complex objectives and (S 2 ^ - j ) whose magnitudes 
are to be minimized. The least square optimization procedure 
discussed in Appendix B is used. The layout for the optimum 
case is shown in Pig. 8.9. Various dimensions of the opti- 
mized circuit in microstrip configuration (f - 1 GHz, 
e r = 10, d = l/4") are also indicated in the figure. 
Effective lengths of 50n and 35.35fi segments are 1.633 cm. 
and 2.243 cm. respectively. In this design also, the micro- 
strip version is obtained by using the technique of scaling 
discussed in Appendix G. 

The variation of input VSWR with frequency is shown in 
Pig. 8.10. Figure 8.11 shows the variations of isolation 
(i.e. |S^[) and of the transfer coefficients to coupled 
ports (i.e. 1^2 i t an ^ |S^j ) with frequency. The input 

VSWR at the center frequency is 1.019 and the difference in 
the magnitudes of coupling coefficients to the two coupled 
ports at the center frequency is about 0.004. 

8.3 DISCUSSION 

In the optimization of input VSWR of power dividers, the 
objective function contains only . This is because the 
symmetry of the circuit ensures that S 21 “ Sjl and when S^ ^ 






VS WR 



FREQUENCY (GHz) 


of optimized hybrid 


ISOLATION 
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tends to zero the power division objective (i.e. equal power 
in the two output ports) is automatically met,- In the case 
of branch-line hybrid, the only symmetry is that the two line 
sections of same impedance have equal lengths. This does not 
ensure equal power division when input VSWR is minimized to 
unity. For this reason it is necessary to have the objective 

*r 

as minimization of and (S 21 - j S^). This would give 
the optimum which tries to minimize input VSWR as well as 
equalise the power division. 

It has been shown in this chapter that the two- 
dimensional approach can be used in the optimization of strip- 
line and microstrip circuits. Improved designs can be obtained 
for circuits whose performance is limited by the presence 
of the discontinuity effects and the higher order modes 
present thereby. 



Chapter Nine 


C ONC HIDING- REMARKS 

This chapter contains summary of the results reported 
in this thesis and some suggestions for further investigations. 

9.1 SUMMARY OR SALIENT RESULTS 

The investigations reported in this thesis can be broadly 
divided into two categories - i) the analytical techniques 
developed for two-dimensional circuits and ii) their appli- 
cations to microwave integrated circuits. 

9.1.1 Analytical Techniques 

Green's functions for some triangular and sectoral shapes 
have been reported. 

For triangular shapes, the method of images [25] has been 
employed to obtain the eigenfunctions of the Helmholtz equation 
which are then used to derive the Green's functions. The tri- 
angular shapes considered are - i) 30° - 60° right-angled, 
ii) equilateral, and iii) right-angled isosceles. Green's 
functions for odd mode half -sections of some of the triangular 
shapes and of a rectangle are al.so reported. These segments 
have mixed boundary conditions and the method of images is 
applied in these cases also. 
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The Green's functions for circular sectors, annular rings 
and annular sectors have been obtained by expanding these in 
series of eigenfunctions [25J. As mutually orthogonal eigen- 
functions are obtainable for only those sectoral angles which 
are submultiples of it, the Green's functions are obtainable 

for only these sectoral shapes. For circular sectors, the 
2 

eigenvalues m are given by the zeros of the first deriva- 
tives of Bessel’s functions. Bor annular rings and annular 
sectors, the transcedental equation (4.22) has to be solved 
for evaluating various eigenvalues. 

The segmentation method for the analysis of 2-d circuits 
has been modified so that the operations are carried out in 
terms of Z-matrices of individual segments. S -matrix for the 
overall circuit is obtained from its Z-matrix. The proposed 
method using Z-matrices is shown to be computationally more 
efficient than the earlier method. 

9.1.2 Applications 

Various modes of resonances for equilateral triangular, 
60°-rhombic and regular hexagonal resonators have been 
studied. These shapes have been analyzed using Green's fun- 
ction for an equilateral triangle. Various modes of resonances 
are found for these three shapes and the voltage distributions 
around the periphery for these modes have been obtained. It is 
seen that the triangular resonators have odd modes of resonances 



in addition to the even modes reported earlier. The 
resonance frequencies predicted have been verified experimen- 
tally for each of the three shapes. 

Stripline and microstrip T-junctions have been shown to 
exhibit lower parasitic reactances when triangular geometries 
are used as compared to the junctions with rectangular geometry 
normally used. Two types of T-junctions with impedance ratios 
1:1:1 and (l/V”2):1:1 have been studied for this purpose. For 
T-junctions with equal impedances, an equilateral triangle 
and for T-junctions with impedance ratio (l/V"2):1:1, a right- 
angled isosceles triangle is proposed for use at the junction. 
Improvements using triangular junctions have been verified 
experimentally for both the cases. 

The 2-d analysis technique has been used for compensation 
of discontinuity reactances and optimization of power dividers 
and branch-line hybrids. Triangular geometries have been used 
at T-junctions. A least square optimization procedure is used 
for obtaining the designs of these circuits. Experimental 
results for one type of power divider circuit in microstrip 
configuration are also reported. 

9.2 SUGGESTIONS FOR FURTHER INVESTIGATIONS 

Green's functions for some other shapes can also be 
obtained by expansion in series of eigenfunctions. Using the 
eigenfunctions for elliptical shapes. Green's functions for 
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ellipses, elliptical sectors and elliptical rings can be 
derived. 

Triangular, rhombic and hexagonal resonators studied 
can also be used as microstrip antennas by determining the 
radiation field around the periphery. Many other resonator 
shapes (like annular rings, sectors, etc.) can also be 
analyzed by the Green's function approach and used to design 
antenna configurations of similar shapes. 

Triangular T-junctions considered in this thesis could 
also be used for improving the performance of other circuits 
which contain any of the two types of T-junctions discussed in 
Chapter Seven. Some examples of these circuits are SPDT 
switches, stub filters, ratrace hybrids, etc. 

The technique used for minimizing the effect of parasitic 
reactances in T-junctions can be extended to junctions with 
arbitrary impedances also. In these junctions an improvement 
in performance is expected by using triangles with the sides 
equal to the effective widths of the three outgoing lines. 
Contour integral method [1 ] can be used for the analysis of 
these junctions. The T-junctions with arbitrary impedances 
are used in several microwave circuits like multi -branch 
directional couplers, filters, impedance matching stubs, etc. 
Improvements in the performance resulting from the use of 
triangular T-junctions in these circuits should be investi- 
gated . 
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Planar circuit approach has been used in discontinuity- 
analysis £l5]. With the Green's functions for triangular and 
sectoral shapes reported in this thesis, more generalized 
discontinuity configurations can now he analyzed. Seme of 
these configurations are shown in Figures 4.1, 4.2 and 7.1. 
Methods for compensation of discontinuity reactances can be 
studied using 2-d approach. 

The 2-d analysis technique can, in general, be used in 
optimization of the design of microstrip and stripline cir- 
cuits. Considerable improvement is expected for the circuits 
whose performances are limited by the junction reactances and 
the presence of higher order modes. A general CAD program can 
be developed for this purpose. 



Appendix A 


INTEGRALS OP EIGENFUNCTIONS IN GREEN’S 
FUNCTION EXPANSIONS 


In this Appendix, normalized integrals of eigenfunctions 
(appearing in Green's function expansion) over the widths of 
coupling ports are given. Normalized integral is obtained by 
dividing the integral of eigenfunction over the width of 


port by the port width. These integrals have been used for 
the characterizations for rectangular, equilateral triangular 
and right-angled isosceles triangular components in Chapters 
Six, Seven and Eight. In the expressions which follow, fun- 
ction sine (x) is given by 



A. 1 RECTANGULAR SEGMENT 

Normalized integrals of the eigenfunction 
cos cos (given in equation (2.30)) are being given 

for the ports along each of the four sides as shown in 
Fig. A. 1(a). These integrals are as given in Table A. 1 . 

Green's function for odd mode half -section of a rectangle 
is given in Sec. 3.5. This function is used to carry out odd 
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Table A, 1 Normalized integrals for ports on a rectangular 
s egment 

locations of ports Normalized Integrals 

Along the side OC 

Along the side OA 

Along the side AB 
Along the side CB 



Along the side OC sin (®jpO sincC^I^-) 

Along the side AB (-l) m sin(^p-) sinc(^^) 


S=1 

(-1) 2 cos(nms) sinc(^~) 


cos (nus) sinc(^^) 

cos (nms) sinc(^~) 

(-l) m cos(mts) sinc(^^) 
(-l) n cos(mits) sinc(^r^) 


Along the side CB 
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mode analysis of power dividers in Sec. 8.1. The eigen- 
functions are given by cos sin where n takes only 

odd values. The normalized integrals for ports along the 
three sides as shown in Pig. A. 1(a) are given in Table A. 2. 

A. 2 EQUILATERAL TRIANGULAR SEGMENT 

The eigenfunctions of equilateral triangular segment 
(Pig. A.1(b)) are Tj(x,y) and T 2 (x,y) given by (3.5) and 
(3.15) respectively. Normalized integrals (denoted by F) 
of T^(x,y) and T 2 (x,y) are being given for ports along each 
of the three sides as shown in Pig. A. 1(b). In the expre- 
ssions which follow, a and {3 are given by and res- 
pectively. The normalized integrals for T^(x,y) are as given 
below. 

For ports along the side AB : 

P = (-1) 2, cos[p(m-n)] sinc[a(m-n)] 

+ (-l) m cos[p(n-S.)3 sinc[a(n-£)] 

+ (-l) n cos[p(A-m)] sinc[a(£-m)] . (A. 2) 

For ports along the side AC : 

P = cosCp(m-n)] sinc[a(m~n)3 + cos[p(n-£)] sinc[a(n-fc)] 

+cos[p(£-m)] sinc[a(s.-m)] . (A. 3) 

For ports along the side BC , the expression for P is same as 
that for ports along the side AC. The normalized integrals for 
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^ 2 (x,y) are as given below. 

For ports along the side AB : 
t 

F = (-1) sin[p(m-n)J sinc[a(m-n)] 

+ (-l) m sin[p(n-£)] sinc[a(n-£)] 

+ (-l) n sin[p(t-m)] sinc[ a( £-m)] . (A.4) 

For ports along the side AC : 

F = -sin[ p(m-n)] sinc[a(m-n)] - sin[p(n-£)l sine[ a(n-fl. )] 

-sin[ P( S'-m)] sinc[a(£-m)] . (A. 5) 

For ports along the side BC , the expression for F is negative 
of that for the ports along AC. 

A. 3 RIGHT-ANGLED ISOSCELES TRIANGULAR SEGMENT 

The eigenfunctions of right-angled isosceles triangle 
(Fig. A, 1(c)) are given by T(x,y) in (3-24). Normalized 
integrals of T(x,y) are being given for ports along each of 
the three sides as shown in Fig. A.l(c). The normalized 
integrals for T(x,y) are as given below. 

For ports along the side OA, 

F = cos(nixs) sinc(^~) + (-l) m+n cos(mn;s) sinc(^~). (A. 6) 
For ports along the side OB, 

F = cos(tft 7 is ) sinc(^vj~) + (-1) m+n ' cos(mxs) sinc(^ 2 ~). (A.7) 
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For ports along the side IB, 

F = (-1 ) m { cost {m+nj-rcs] sinc[ (ntf-n)aj + 

cos[ (m-n)ns] sinc[(m-n)a] > (A.8) 

where 

a = . 

2f5a 

Green’s function when side OB in Fig. 1.1(e) is an 
electric wall is given by (3.29). This function is used to 
carry out odd mode analysis of power dividers in Sec. 8.1. 

The eigenfunction U(x,y) is given by (3.30). Normalized 
integrals, F, of U(x,y) for ports along two sides of the 
triangle are given below. 

For ports along the side OA, 

m+n 

F = sin(*vp) sinc( s ||-) - (-1) 2 sin( 2 | s ) sint 2 ^) . (A.9) 

For ports along side AB, 
hh- 1 

F = (-1) 2 {cos[ (cH-n)p] sinc[ (mfn)a] - 

cost (m-n)p] sinct (m-n)a] > (A. 10) 

where 

a ~ and |3 = ns/2 . 

4Y2a 



Appendix B 


OPTIMIZATION PROCEDURE USED 

This appendix describes the procedure used for optimization 
of circuits in Chapter Eight. The method uses the gradient of 
the performance objectives which are expressed in terms of S- 
parameters. Differential Z-matrices of components with respect 
to the designable parameters are obtained by series summation 
in the same manner as used for Z-matrices of components. Por 
segmentation, these differential Z-matrices are put together 
in the same form as in equation (5.11). The differential Z- 
matrix of the overall network is obtained using (5.14) as 

K - V + - V - + W 1 <4p - V 

" ^pr^^qq ~ \r * ^rq + ^^qq ~ ^qr ~ ^rq + - rr ^ 

• (Z, q - V - * Sr*.)'’ <2qp - %> + <»•«> 

The differential Y-matrix is obtained using 

11 = “C 1 K si 1 ( B ‘2) 

-p ~P ~P -p 

and the external subports are combined in the same manner as 
in (5.23). The differential S-matrix of the overall network 
is obtained using (5.24) as 

s' = -[fz 0 Ki+ a, - 1) dp + x)*’] i' (I, + X)-' [mg.(B.3) 
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Having obtained the differential S-matrices, the 
gradient of performance can be determined and any of the 
gradient optimization methods can be used. A least square 
optimization procedure is utilized in Chapters Seven and 
Bight. 

Optimization of least Square Ob jective ^F unctions 

The procedure described here minimizes the sum of squares 
of deviations. It is based on Gaussian least square solution. 
The method uses the gradients of each of the terms and avoids 
one -dimensional search in the design space. 

Let the objective function U to be minimized be written 
as 

U = f^ + f\ + ... + f^ (B.4) 

where f ^ ...» f m are the functions of the designable 
parameters x given by 

£. = { x ^ » %2 > * t ( B » 5 ) 

with m>_n. The minimization of U in (B.4) is equivalent to 
minimizing the norm of an m-dimensional objective function 
vector 

f(x) ='(f 1 (x) > f 2 (x), ...» f m (x) f (B.6) 

with tA&), f 2 (s), being its components. The norm 

of f(^) is taken to be Cf(x)]^ [f(x)] and is equal to TJ given 
by (B.4). The Jacobian, formed using derivatives of the 
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components of the vector objective function, is given as 



(B.7) 


The value of the vector objective function f can be est ima ted 
at x + ax using the Taylor series expansion as 


+ A &) * £(s) + l(x) ax 


(B.8) 


where the increment Ax is assumed to be small so that the 
higher order terms can be neglected. 

In this method, Ax is chosen such that the norm of 
£.(& + A 2L)> given by the right-hand side of (B.8), is minimum. 

The value of Ax giving minimum norm of f(x) + J(^) Ax is 
given by [34-] 

*S = - (Cite)]* [its )])' 1 tits)]* its) . (B.9) 

The process converges in a single iteration if U given by (B.4)I 
quadratic. Bor higher order objective functions, this process 
has to be repeated i.e. £(&) and J(£) are computed at x + Ax 
and a new A& is obtained. This process is repeated till 


convergence. 
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A common drawback using this method is that A£, as 
computed from (B.9) may become too large. In such cases it 
is desirable* to restrict each component of as to a pre- 
specified value. 

This technique has been used for optimization in this 
thesis. The scattering parameters, being complex, are treated 
as two components of the vector objective function. For 
example, in the case of minimizing input reflection coefficient, 
f is expressed as 

f = {Re(S t1 ), ImCS^))*. (B.10) 

In the design of branch-line hybrid the objectives for optimi- 
zation are to minimize the input VSWR and to equalize power 
division. Objective function vector f is expressed as 

£ = {Re(S n ), Im(S n ), Re(S 21 -j S 41 ), Im(S 2l - j . 

(B.11) 

The results reported in Chapters Seven and Bight are based on 
the optimization procedure described above. 



Appendix C 


SCALING FOR PLANAR CIRCUITS 


Once a 2-d circuit has been designed for a given set of 
specifications, the principle of scaling may be used to design 
a circuit with similar specifications at a different frequency 
or with different (scaled) impedances at all the ports- Design 
of a stripline type circuit from a microstrip type circuit (and 
vice versa) is also possible through scaling. 

C.1 EFFECTIVE DIMENSIONS 


For the purpose of scaling, the dimensions considered 
include the extension of periphery for taking into account 
the fringing field in open boundary planar circuits. In 
stripline type 2-d circuits, the effective width of a trans- 
mission line of impedance Z Q is given by 


W 


eff 



(C.1) 


and equals (W + 26) with 

c 2d *n 2 

In microstrip type 2-d circuits, the effective width of 
transmission line of impedance Z 0 is 


y — XUl 
w eff Z Q * 

The width W is related to the impedance Z Q by [35] 


(C.3) 



W = 8d 


{ C exp ( Z Ve"+?/ 42 * 4 ) -1 J ( 7+4/ e )/ 1 1 + (1+l/e J/0.81 }‘ 


o' r 


exp ( Z Q V / 42 . 4 ) - 1 


and 6 is computed as 

<5 = (W eff - W)/2 


(C.4) 


(C.5) 


Consider a rectangular segment with effective length and 

width £ and W respectively. In strip line configuration 6 
0 0 

is almost independent of width and the extensions along the 

length and width are same. For microstrip configuration, 

first Zo is computed as r\ d/ s, . The actual length £ is 

0 

computed from Z^ using (C.4). The extension along the length 
is obtained as ( a -a )/2. The extension along the width is 
computed as described above using (C.3) to (C.5). 

C .2 FREQUENCY SCALING 

The concept of frequency scaling applicable to strip line, 
microstrip and waveguide circuits [36], [37] can be applied 
to general 2~d planar circuits also. If all the linear dimen- 
sions in the planar patch are reduced by a scaling factor E 
and it is deposited on a substrate of height d/K, the response 
of the modified circuit at frequency f will be same as that of 
the original circuit on a substrate of height d at frequency 
f/K. For example, if a circuit has been designed at 1 GHz on 
a 5 mm thick substrate, by scaling down all the linear dimen- 
sions of the patch by a factor of 10 and fabricating It on a 
0.5 mm thick substrate the design can be used at 10 GHz, 
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It may be noted that for strip line type circuits, exten- 
sion of boundary because of the fringing field is proportional 
to the height of the substrate. Thus, actual physical dimen- 
sions may be used for the purposes of frequency, sca ling - 

C.3 IMPEDANCE SCALING 

The concept of impedance scaling gives another planar 
circuit which has impedance values at all the coupling ports 
scaled by some factor and has the same frequency response for 
the scattering parameters. The new circuit has to be deposited 
on a substrate of different height and all dimensions of the 
2-d patch remain the same. If the planar patch is deposited 
on a substrate of height Kd, the S-matrix of the new circuit 
would be the same as that of the original circuit on substrate 
height d. For the new circuit, the impedances at all the 
coupling ports are scaled up by the factor K. 

The two types of scaling concepts discussed above may be 
combined together. For example, if the height of the substrate 
is kept unchanged and linear dimensions scaled up by a factor K, 
the S-matrix of the modified circuit at the frequency f/K is 
the same as that of the original circuit at frequency f but all 
impedances levels (in the modified circuit) are scaled down by 
the factor K. 

C.4 DESIGN OF A STRIPLINE TYPE CIRCUIT FROM A MICROSTRIP TYPE 

PLANAR CIRCUIT (AND VICE VERSA) 

The concept of scaling can also be applied to obtain an 
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equivalent design in the micros trip type configuration from the 
design in the stripline type configuration and vice versa. The 
scattering matrix response of a etripline type planar circuit 
with spacing 2d between the ground planes will be the same as 
that of a microstrip type planar circuit with substrate height 
d, provided the linear dimensions (between the magnetic walls) 
in the 2-d patches are the same in both the cases. It may be 
be noted that the characteristic impedance values at the coupling 
ports for the microstrip type circuit will be twice the corres- 
ponding values for the stripline type circuit. Impedance scaling 
technique can be used in conjunction with the above conversion 
to get the desired impedance values at the coupling ports. 

These scaling techniques provide an additional flexibility 
in design and are used to obtain the designs of circuits in 
microstrip configuration in Chapter Eight, 
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